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ABSTRACT 

We study the shapes of galaxy dark matter haloes by measuring the anisotropy of the weak gravitational lensing signal 
around galaxies in the second Red-sequence Cluster Survey (RCS2). We determine the average shear anisotropy within 
the virial radius for three lens samples: the 'all' sample, which contains all galaxies with 19 < m r i < 21.5, and the 
'red' and 'blue' samples, whose lensing signals are dominated by massive low-redshift early-type and late-type galaxies, 
respectively. To study the environmental dependence of the lensing signal, we separate each lens sample into an isolated 
and clustered part and analyse them separately. We address the impact of several complications on the halo ellipticity 
measurement, including PSF residual systematics in the shape catalogues, multiple deflections, and the clustering of 
lenses. We estimate that the impact of these is small for our lens selections. Furthermore, we measure the azimuthal 
dependence of the distribution of physically associated galaxies around the lens samples. We find that these satellites 
preferentially reside near the major axis of the lenses, and constrain the angle between the major axis of the lens and 
the average location of the satellites to (6) = 43.7° ± 0.3° for the 'all' lenses, (6) = 41.7° ± 0.5° for the 'red' lenses 
and (6) = 42.0° ± 1.4° for the 'blue' lenses. We do not detect a significant shear anisotropy for the average 'red' and 
'blue' lenses, although for the most elliptical 'red' and 'blue' galaxies it is marginally positive and negative, respectively. 
For the 'all' sample, we find that the anisotropy of the galaxy-mass cross-correlation function (/ — /4s) = 0.23 ± 0.12, 
providing weak support for the view that the average galaxy is embedded in, and preferentially aligned with, a triaxial 
dark matter halo. Assuming an elliptical Navarro-Frenk- White (NFW) profile, we find that the ratio of the dark matter 
halo ellipticity and the galaxy ellipticity fh = en/e g = 1.50i}'o 3 , which for a mean lens ellipticity of 0.25 corresponds to 
a projected halo ellipticity of = 0.381q ,2 ,5 if the halo and the lens are perfectly aligned. For isolated galaxies of the 'all' 
sample, the average shear anisotropy increases to (/ — /4s) = 0. 51 jlo'lf and fh = 4.731?, 05, whilst for clustered galaxies 
the signal is consistent with zero. These constraints provide lower limits on the average dark matter halo ellipticity, as 
scatter in the relative position angle between the galaxies and the dark matter haloes is expected to reduce the shear 
anisotropy by a factor ~2. 

Key words, gravitational lensing - dark matter haloes 



1. Introduction 

Over the last few decades a coherent cosmological 
paradigm has developed, ACDM, which provides a frame- 
work for the study of the formation and evolution of 
structure in the Universe. N-body simulations that are 
based on ACDM predict that (dark) matter haloes collapse 
such that their density pro files closely follow a Navarro- 
Frenk- White profile (NFW; iNavarro et al.lfl996ri . which is 
in excellent agreement with observations. Another fun- 
damental predic tion from simulations is that the haloes 
are t riaxial (e.g. iDubinski fc Carlberdll99ll lAllgood et al.l 
[200l . which appear elliptical in projection. This predic- 
tion of dark matter haloes, as well as many others con- 



cerning the evolution of their shapes (e.g IVera-Ciro et al.l 
1201 it ) . the effect of the central galaxy on the dark mat- 
ter halo shape (e.g. Kazantzidis et al.l l2010t lAbadi et al.l 
120101: iMachado fc Athanassoulal I2010D and the ir depen- 



dence on environment (e.g Wang et al.l 1 2 1 lh . remain 
largely untested observationally. 

Direct observational constraints on the halo ellipticitics 
have proven to be difficult, mainly due to the lack of useful 
tracers of the gravitational potential. On small scales (~few 
kpc) , halo ellipticity estimates have been obtained through 
the combination of strong lensing and stellar dynamics 
fe.g.lvan de Ven et al.lkoialDutton et al.ll201lHSuvu et al.l 
12011 ). planetary nebulae (e.g. iNapolitano et al.ll201 If) and 
HI observations in late-type galaxies fe.g. iBaneriee fc Jog] 
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120081: lO'Brien et alJ[2010h . On larger scales, the distribu- 
tion of satellite ga laxies around centrals has been used (e.g. 
iBailin et al.ll2008l ) , but such studies have only provided con- 
straints for rich systems that may not be representative for 
the typical galaxy in the universe. 

Weak gravitational lensing does not depend on the pres- 
ence of optical tracers and is capable of providing ellipticity 
estimates on a large range of scales (between a few kpc to 
a few Mpc). Therefore it is a powerful observational tech- 
nique to study the ellipticity of dark matter haloes. In weak 
lensing the distortion of the images of faint background 
galaxies due to the dark matter potentials of intervening 
structures, the lenses, is me asured. This has been used to 
determine halo masses (e.g. Ivan Uitert et al.l [20111) as well 
as the extent of haloes. If galaxies preferentially align (or 
anti-align) with respect to the dark matter haloes in which 
they are embedded, the lensing signal becomes anisotropic. 
This signature can be used to c onstrain the ellipticity of 
dark matter haloes of galax ies (jBrainerd fe Wright] l2000t 
iNataraian fc Refregierl 1200 Oft . 

The core assumption in the weak-lensing-based halo 
ellipticity studies is that the orientation of galaxies and 
dark matter haloes are correlated; if they are not, the 
shear signal is isotropic and cannot be used to constrain 
the ellipticity of the haloes. The relative alignment be- 
tween the baryons and the dark matter has been addressed 
in a large number of studies based on numerical simu- 
lations fe.g. Ivan den Bosch et al.l l200l l2003t IBailin et all 



20051: iKang et all 120071: iBett" et al. M2010l: lHahn et al.ll2010t 



Deason et al. 20111) . in studies based on the distribution 



of satellite galaxies around centrals (jWang et al.l l2008t 
lAgustsson fc Brainerdl [2010h an d in studies based on the 
ellipticity correlation function ([Faltenbacher et ahl 120091 : 
lOkumura et al.l l2009t) . The general consensus is that al- 
though the galaxy and dark matter are aligned on average, 
the sc atter i n the differential position angle distribution is 
large. iBettl (|201ll) examined a broad range of galaxy-halo 
alignment models by combining TV-body simulations with 
semi-analytic galaxy formation models, and found that for 
most of the models under consideration, the stacked pro- 
jected axis ratio becomes close to unity. Consequently, the 
ellipticity of dark matter haloes may be difficult to measure 
with weak lensing in practice. 

Knowledge of the relative alignment distribution is 
not only crucial for halo ellipticity studies, but also for 
studies of the intrinsic alignments of galaxies. Numerical 
simulations predict that the shapes of neighbouring 
dark matter haloes are correlated (e.g . ISplinter et al 



1997 tlCroft fc Metzlerll2000l: iHeavens et al.ll2000l: iLee et al 



20081) . The shapes of galaxies that form inside these haloes 



may therefore be intrinsically aligned as well. Measuring 
this effect is interesting as it provides constraints on struc- 
ture formation. Also, the lensing properties of the large- 
scale structure in the universe, known as cosmic shear, are 
affected by intrinsic alignments, and benefit from a careful 
characterization of the effect. Intrinsic alignments arc stud- 
ied observationally by correlating the cllipticities of galax- 
ies as a function of separation; misalignments can signif- 
icantly reduce these ellipticity correlation functions (e.g. 
iHevmans et"al"1l2004l) . 

To date, only three observational weak lensing 
studies have detected the anisotropy of the lensing 
signal dHoekstra et all 12004 iMandelbaum et al.l l2006at 
iParker et al.ll2007l) . These studies have provided only ten- 



tative support for the existence of elliptical dark mat- 
ter haloes, as they were limited by ei ther their sur- 
vey size and lack of c olour informati on dH oekstra et al 



2004 IParker et al.ll2007l) or their depth (jMandelbaum et al 
2006afl . To improve on these constra ints, we use the Red 



sequence Cluster Survey 2 (RCS2, iGilbank et aD[20ll . 
Covering 900 square degree in the g'r' z' -bands, a limit- 



ing magnitude of 



lim 



24.3 and a median seeing of 



0.7", this survey is ve ry well suited for lensing studies (see 
Ivan Uitert et al.ll201hY Using the colours we select massive 
luminous foreground galaxies at low redshifts. To investi- 
gate whether the formation histories and environment affect 
the average halo ellipticity of galaxies, the lenses are sepa- 
rated by galaxy type and environment, and the signals are 
studied separately. 

The structure of this paper is as follows. We describe the 
lensing analysis, including the data reduction of the RCS2 
survey, the lens selection and the definition of the shear 
anisotropy estimators, in Section [5] We present measure- 
ments using a simple shear anisotropy estimator in Section 
and use it to study the potential impact of PSF resid- 
ual systematics in the shape catalogues. Various compli- 
cations exist that might have altered the observed shear 
anisotropy, and in Section 0] we study the impact of two of 
them: multiple deflection and the clustering of the lenses. 
The shear anisotropy measurements are shown and inter- 
preted in Section [5] We conclude in Sectio nal Throughout 
the p aper we assume a WMAP7 cosmology (jKomatsu et al.l 
[20Tlh with cr 8 = 0.8, Ov = 0.73, fl M = 0.27, Sl b = 0.046 
and the dimensionless Hubble parameter h = 0.7. The er- 
rors on the measured and derived quantities in this work 
generally show the 68% confidence interval, unless explic- 
itly stated otherwise. 



2. Lensing analysis 

For our lensing analysis we use the imaging data 
from the second R ed-sequence Cluster Survey (RCS2; 
IGilbank et alll201lh . The RCS2 is a nearly 900 square de- 
gree imaging survey in three bands (</ , r 1 and z 1 ) carried out 
with the Canada- France-Hawaii Telescope (CFHT) using 
the 1 square degree camera McgaCam. In this work, we use 
the ~700 square degrees of the primary survey area. The 
remainder constitutes the 'Wide' component of the CFHT 
Legacy Survey (CFHTLS) which we do not consider here. 
We perform the lensing analysis on the 8 minute exposures 
of the /-band {r' lim ~24.3), which is best suited for lensing 
with a median seeing of 0.71". 

2.1. Data reduction 

The phot ometric calibra t ion of the RCS2 is described 
in detail in IGilbank et al.l (|2011l) . The magnitudes are 
calibrated using the colours of the stellar locus and the 
overlapping Two-Micron All-Sky Survey (2MASS), and 
have an accuracy better than 0.03 mag in each band 
compared to the SDSS. The crea tion of the galaxy shap e 
catalogues is described in detail in Ivan Uitert et aU(|201ll ). 
We refer readers to that paper for more detail, and present 
here a short summary of the most important steps. 

We retrieve the ElixiiQ processed images from the 



http:/ /www. cfht.hawaii.edu/Instruments/Elixir/ 
2 http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 
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Canadian Astronomy Da ta Centre (CADC ) arch ivfl Wc 
use the THELI pipeline (|Erben et al.ll2005l l2009h to sub- 
tract the image backgrounds, create weight maps that we 
use in the object detection phase, and to identify satellite 
and asteroid trail s. To detect the objects i n the images, we 
use SExtractor (jBertin fc Arnoutslll99ll . The stars that 
are used to model the PSF variation across the image are 
selected using size-magnitude diagrams. All objects larger 
than 1.2 times the local size of the PSF arc identified as 
galaxies. We measure the shapes of the galaxies with the 
KSB method ([Kaiser et all 119951: iLuppino fc Kaiserlll997t 
IHoekstra et al. 19981 ) . usin g the im plementation described 
bv IHoekstra et all (|1998L |2000() . This implementation 
has been tested on simulated images as part of the 
Shear Testing Programmes (STEP) (the 'HH' method in 
Hcymans et al. 2006 and Massey et al. 2007), and these 
tests have shown that it reliably measures the unconvolved 
shapes of galaxies for a variety of PSFs. Finally, the 
source ellipticities are corrected for camera shear, which 
originates from slight non-linearities in the camera optics. 
The resulting shape catalogue of the RCS2 contains the 
ellipticities of 2.2xl0 7 galaxies. A more detailed discus - 
sion of the analysis can be found in Ivan Uitert et all (|2011l) . 



2.2. Lenses 

To study the halo ellipticity of galaxies, we measure the 
shear anisotropy of three lens samples. The first sample 
contains all galaxies with 19 < m r > < 21.5, and is referred 
to as the 'all' sample. This sample consists of different 
types of galaxies that cover a broad range in luminosity 
and redshift. The shear anisotropy measurement of this 
sample enables us to determine whether galaxies are on 
average aligned with their dark matter haloes. 

The formation history of galaxies differs between galaxy 
types, and consequently the relation between baryons and 
dark matter may differ too. Therefore, the average dark 
matter halo shapes, and the orientation of galaxies within 
these haloes, might depend on galaxy type. To examine 
this, we separate the lenses as a function of their type. 

Various selection criteria have been employed 
to separate early-type from late-type galaxies. In 
most cases, galaxies are either sele cted based on the 
slope o f their brightne s s pro files (|Mandelbaum et all 



ia; 



i 2on. 



or on 



Hoekstra et al. 



the ir colours 
20051). To study 



I2006bl Ivan Uitert et al 
(jMandelbaum et al.ll2006 

how th ese selection criteria relate, Mandelbaum et ahl 
(l2006aft compare the selection based on their SDSS u — r 
model colour to the selection based on the frac-dev 
parameteiB, and find that the assigned galaxy types agree 
for 90% of the galaxies. 

We choose to separate the galaxy types based on 
their colours, as the g'-, r'- and z'- band colours arc 
readily available for all galaxies in the RCS2. The aim 
of the separation is two-fold: to make a clean separation 
between the red quiescent galaxies which typically exhibit 
early-type morphologies and blue star-forming galaxies 
that typically have late-type morphologies, and to select 



The jracjdev parameter is determined by simultaneously fit- 
ting frac-deV times the best-fitting de Vaucouleur profile plus 
(1-fracjdeV) times the best-fitting exponential profile to an ob- 
ject's brightness profile 



massive lenses at low redshifts to optimize the lensing 
signal-to-noise, and minimize potential contributions 
from multiple deflections (see Section 14. ip . To determine 
where these massive low-redshift galaxies reside in the 
colour-magnitude plane, we use the photometric redshift 
catalogues of the CFHTLS Wide from the CFHTLenS 
collaboration (jHildebrandt et all 1201 lh . and define our 
boxes accordingly; details of the selection of the 'red' and 
'blue' lens sample are described in Appendix [A] Note that 
these lens samples overlap with the 'all' sample, but not 
with each other. Details of the samples are given in Table 

m 

To study how well we can separate early-types from 
late-types, we compare our selection to previously em- 
ployed separation criteria. Details of the comparison can be 
found in Appendix [A] Wc find that the 'red' sample is very 
similar to the selection based on the u' — r' colour, whilst 
^58% of the 'blue' sample are actually red according to 
their u' — r' colour. Most of these contaminants of the 'blue' 
sample are not massive, and actually dilute the lensing 
signal. The purity of the 'blue' sample could be improved 
by shifting the selection boxes to bluer colours, but this at 
the expense of removing the majority of massive late- type 
lenses. Finally, we note that ^70% of the 'all' sample are 



considered blue based on their 



colours. 



To study the second objective of the lens selection, i.e. 
to select massive and bright low-redshift lenses, wc apply 
the colour cuts to the CFHTLS Wl photometric catalogue, 
and show the distribution of absolute magnitudes and 
photometric redshifts of the lens samples in Figure [T] We 
find that the 'red' lens sample consists of galaxies with 
absolute magnitudes in the range — 24 < M T < —22, and 
most with redshifts between 0.3 and 0.6. The galaxies 
from the 'blue' sample have absolute magnitudes in the 
range —24 < M r < —20, and are located at redshifts 
between 0.1 and 0.6. For the blue galaxies, we cannot 
define a criterion that exclusively selects luminous lenses 
in a narrow redshift range, based on the g'r'z' magnitudes 
alone. Finally, the 'all' sample has the broadest luminosity 
and redshift distribution. It is possible to narrow down 
the redshift range by discarding the lenses with the largest 
apparent magnitudes from each sample. We choose not 
to, however, because this lowers the signal-to-noise of the 
lensing measurement, which consequently broadens the 
constraints on the average halo ellipticity. 

Note that due to the lack of a very blue observing band 
in the CFH TLS, the photometric red shifts below 0.2 are 
biased high (jHildebrandt et al"1l201l|) . As a consequence, 
a fraction of the galaxies of the 'blue' lens sample may 
have been shifted to higher redshifts, and thus larger 
luminosities. The mean redshift and luminosity of the 
sample may therefore be somewhat smaller than the values 
quoted in Table [TJ and the distributions shown in Figure [T] 
are only indicative. 

Since the dark matter halo ellipticity is measured 
relative to the ellipticity of the galaxy, it is interesting 
to examine the distribution of the latter. In Figure [2J we 
show the ellipticity distribution of the lens samples; the 
mean galaxy ellipticity of each sample is given in Table [TJ 
The ellipticity distributions of the 'all' and 'blue' sample 
are comparable, and are broader than the 'red' sample 
one, because the 'all' and 'blue' sample have a considerable 
fraction of disk galaxies. The differences between the el- 
lipticity distributions have consequences for the weighting 
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M , Photometric redshift 

r 

Fig. 1. Number of lenses as a function of absolute magnitude (a) and redshift (b) for the three lens samples, ob- 
tained by applying ident ical cuts to the CFHTLS Wl photometric redshift catalogue from the CFHTLenS collaboration 
(jHildebrandt et al1l201l|) . The 'all' sample (black solid lines) has the broadest distributions, and covers absolute r'-band 
magnitudes between —18 and —24, and redshifts between and 0.6. The luminosities of the 'blue' sample (blue dotted 
lines) are in the range —24 < M r < —20, with redshifts 0.15 < z < 0.6. The 'red' sample (purple dot-dashed lines) has 
the narrowest distributions, with luminosities —24 < M r < —22 and redshifts 0.3 < z < 0.6. 

Table 1. Properties of the lens samples: the number of lenses, the mean redshift, the mean luminosity, the mean ellipticity, 
the fraction of lenses that are isolated, the virial mass, the virial radius and the scale radius. 



Sample 


Nl CIls (z) 


(V) 
[10 10 h n 2 L Q ] 




/iso 


M200 
[10 10 h^M e ] 


T200 

[/i^kpc] 


r a 

[/i^kpc] 


All 

Red 
Blue 


1681826 0.31 
136196 0.43 
147079 0.31 


2.86 
8.91 
4.68 


0.25 
0.20 
0.26 


0.20 
0.41 
0.55 


2i.ii?j 

138 ± 8 
44.4±g;| 


150^3 
280 ±5 
192 ±5 


23.9t°; 2 5 
54.6±i;S 

00 7-+0.8 
oz " ' -0.9 



scheme of the lensing anisotropy measurements, as we 
will discuss in Section 12.31 In the analysis, we only use 
galaxies with 0.05 < e g < 0.8, which excludes round lenses 
that do not have a well-defined position angle, and very 
elliptical galaxies whose shapes are potentially affected by 
neighbours and/or cosmic rays. 

The ellipticity of dark matter haloes may depend on 
the environment of a galaxy. We therefore divide the 
lens samples further into isolated and clustered ones, and 
study the lensing signal separately. As we lack redshifts 
for all the galaxies, we have to use an isolation criterion 
based on projected angular separations: if the lens has a 
neighbouring galaxy within a fixed projected separation 
that is brighter (in apparent magnitude) than the lens, 
it is selected for the clustered sample. If the lens is the 
brightest object, it is selected for the isolated sample. We 
test various values for the fixed minimum separation, and 
compare the tangential shear at large scales in Appendix 
[Bl Based on these results, we use a minimum separation 
of 1 arcmin. Note that an environment selection based 



on apparent magnitudes cannot be very pure; a fraction 
of the lenses from the isolated sample may still be the 
brightest galaxy in a cluster. Some of the lenses of the 
clustered sample may in reality be isolated, but have 
been selected for the clustered sample due to the presence 
of bright foreground galaxies. However, the difference 
between the large-scale lensing signal of the isolated and 
the clustered sample indicates that our selection criterion 
works reasonably well. The fraction of the lens sample that 
is isolated, /i so , is indicated in Table [T] 

2.3. Shear anisotropy 

The lensing signal is quantified by the tangential shear, 
7t, around the lenses as a function of projected separation. 
As the distortions are small compared to the shape noise, 
the tangential shear needs to be azimuthally averaged over 
a large number of lens-source pairs: 

<<*>(r) = ^, (1) 

^crit 
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Fig. 2. Ellipticity distribution of the g'r'z'-colour selected 
lens samples. The dashed lines indicate the ellipticity cuts 
we apply to exclude the roundest and most elliptical lenses. 
The ellipticity distributions of the 'all' and the 'blue' sample 
are similar, but the 'red' sample contains relatively more 
round galaxies. 



where AE(r) = £(< r) — S(r) is the difference between 
the mean projected surface density enclosed by r and the 
mean projected surface density at a radius r, and S cr ; t is 
the critical surface density: 



J crit 



4ttG DiDu ' 



(2) 



with Di, D s and Di s the angular diameter distance to 
the lens, the source, and between the lens and the source 
respectively. Since we lack redshifts, we select galaxies 
with 22 < m r i < 24 and a reliable shape estimate as 
sources. We obtain the approximate source redshift distri- 
bution by applying identical magnitude cuts to the pho- 
tometric redshift catalogues of the Canada- France-Hawaii- 
Telescope Legacy Survey (CFHTLS) "Deep Survey" fields 
(jllbert et al.ll2006 l). and find a median source redshift of 
z a = 0.74. This redshift distribution is not exactly identical 
to the one of the sources due to the additional shape param- 
eter cuts applied to the source sample, which are weakly 
dependent on apparent magnitude, but the difference is 
negligble. To convert the tangential shear to AE, we use 
the average critical surface density that is determined by 
integrating over the source redshift distribution: 



1 



J crit/ 



4ttG A 



dz s p(z s ) 



D, 



norm J z\ 
A 



DiD, 



(3) 



dz s p(z s ), 



with p(z s ) the redshift distribution of the sources, and Zi 
the mean redshift of the lens sample used to determine 




Fig. 3. Schematic of a lens galaxy. The tangential shear is 
measured in regions A and B, the cross shear is measured 
in regions C and D. The cross shear is subtracted from the 
tangential shear to correct for systematic contributions to 
the shear. 

Di and Di s . We also measure the cross shear, j x , the 
component of the shear in the direction of 45° from the 
lens-source separation vector. The azimuthally averaged 
cross shear signal should vanish since gravitational lensing 
does not produce it. If this signal is non-zero, however, 
it indicates the presence of systematics in the shape 
catalogues. As the lenses are large and their light may 
contaminate the lensing signal near the lenses, we only 
consider the signal on scales larger than 0.1 arcmin for 
lenses with m r ' > 19, and scales larger than 0.2 arcmin 
for lenses with m r i < 19. These criteria are based on the 
reduction of the source number density near the lenses, 
as discussed in Appendix [D] Hence the smallest scales we 
probe is 28 kpc for the 'all' and 'blue' sample, and 34 kpc 
for the 'red' sample at the mean lens redshift. To remove 
contributions of systematic shear (from, e.g., the image 
masks), we subtract the signal computed around random 
poin ts from the signal co mputed around the real lenses 
f see Ivan Uitert et al.l[20TTb . 

The lensing signal around triaxial dark matter haloes 
has an azimuthal dependence. If galaxies are preferentially 
aligned or oriented at a 90° angle (anti-aligned) with re- 
spect to the dark matter distribution, the lensing signal 
along the galaxies' major axis is respectively larger or 
smaller than along the minor axis, and this dependence 
can be determined. 

To measure the anisotr opy in the si g nal, w e first fol- 
low the approach used by iParker et al.l (|2007|) . For each 
lens, the tangential shear is measured separately using the 
sources that lie within 45° of the semi-major axis (7*3), 
and using those that lie within 45° of the semi-minor axis 
(7t, a) (indicated by B and A in Figure OS respectively). The 
ratio of the shears captures the anisotropy of the signal: 



fmm(r) 



7t,B( r ) 

7*,aW 



(4) 



A value of / mm that is significantly larger (smaller) than 
unity at small scales indicates that the dark matter haloes 
are (anti-)aligned with the galaxies. Systematic contribu- 
tions to the shear, however, may bias the anisotropy of the 
lensing signal. If the systematic shear is fairly constant on 
the scale s where we measure the sig nal, it can be removed 
following iMandelbaum et al.l (|200 6af). In this approach, the 
cross shear component computed in the regions that are 
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rotated by 45° with respect to the major/minor axes (re- 
gion C and D in Figure [3]), 7x,c-d = (7x,c - 7x,d)/2, is 
subtracted from the tangential shear. Spurious shear sig- 
nals contribute equally to jt,A-> 7t,B and 7 x ,c-d, and are 
therefore removed. The corrected ratio then becomes: 



f co "(r) 

J mm V / 



7t,B( r ) +lx,C-p(r) 
7t,A( r ) ~ 7x,C-D(r) ' 



(5) 



If 7x.c-D is zero, the errors on /™m approximately increase 

by a factor v/l + \j\f2\ if 7x.c-d is non-zero, however, the 
errors of f£g£ can either become larger or smaller than 
those of / mm . 

Alternatively, we can assume that the differential sur- 
face density distribution can be descri bed by an isotropic 
part p lus an azimuthally varying part (|Mandelbaum et aLl 
l2006al) : 

(r) = A£ iso (r ) [1 + 2fe g cos(2 AO)} , (6) 

where e g is the observed ellipticity of the lens, A9 is the an- 
gle from the major axis, and / is the ratio of the amplitude 
of the anisotropy of the lensing signal and the ellipticity of 
the galaxy, which i s the p arameter we want to determine. 
IMandelbaum et all (|2006al ) show that the azimuthally vary- 
ing part is given by: 



/A£ iso (r) 



J2i w l AT ll e g ^ i cos(2A6>i) 
2E i ^e 2 „ 4 cos 2 (2A^) ; 



(7) 



with i the index of the lens-source pairs, Wi the weight ap- 
plied to the ellipticity estimate of each source galaxy, which 
is calculated from the shape noise, and e g ^ the ellipticity of 
the lens. This ellipticity is also determined using the KSB 
method, and it is a measure of (1 — i? 2 )/(l + R 2 ) with R 
the axis ratio (R < 1) if the lens has elliptical isophotes. To 
remove contributions from systematic shear, we also mea- 
sure 



/45ASi, 



J2i w t AEi A5 e g ^ cos(2A6>j + ir/2) 
2£ i ^e 2 4 cos 2 (2A^ + 7 r/2) : 



(8) 



where Ei.45 is the projected surface density measured 
by rotating the source galaxies by 45°. The systematic 
shear corrected halo ellipticity estimator is then given by 
(/ - /4 5 )AS iso (r). The average values of / mm , and 
(/ — ^45) within a certain range of projected separations 
arc determined by calculating the ratio of two measure- 
ments for each radial bin, and subsequently averaging that 
ratio within the range of interest. We assume that the er- 
rors of each measurement are Gaussian. Consequently, the 
probability distribution of the ratio is asymmetric, which 
we have to account for. We describe how to calculate the 
mean and the errors of the ratio for a radial bin, and how 
to average that ratio within a certain range of projected 
separations, in Appendix [C] Note that to convert /, the 
anisotropy in the shear field, to //, = eh/e g , the ratio of the 
ellipticity of the dark matter halo and the ellipticity of the 
galaxy, we have to adopt a density profile (e.g. // fh=0-25 
for a singular isothermal ellipsoid, see Mandelbaum et al. 
2006a). 

It is clear from Figure [2] that the ellipticity distribu- 
tions of the red and blue lens samples arc different. It is 
unclear, however, whether the underlying ellipticity distri- 
bution of the dark matter haloes differs as well. If the un- 
derlying distribution is similar for both samples, the pro- 
jected dark matter halo ellipticity cannot depend linearly 



on the galaxy ellipticity. Hence Equation (O might not be 
optimal, and could depend differently on e g . We therefore 
generalise Equation ([7]) to 



/A£ iso (r) = A 



A 



WjAZjeli cos(2Afl t ) 
2£ l ^e 2 «cos 2 (2A# 4 ) 



(9) 



(10) 



and calculate it for different values of a. Equation ([8]l 
changes similarly. The factor A in Equation © scales each 
measurement of .f to the 'standard' of a = 1 as used in 
IMandelbaum et al.l (|2006al) . which cases a comparison of / 
for different values of a. The optimal weight results in the 
best signal-to-noise of the measurement. 

The different halo ellipticity estimators can in princi- 
ple be used to study the relation between the ellipticity of 
the galaxy and the ellipticity of their dark matter hosts. 
In particular, Equation ([5]) is defined such that it depends 
on the average dark matter halo ellipticity, whilst Equation 
Q is sensitive to the relation between the galaxy elliptic- 
ity and the dark matter ellipticity. Hence by comparing the 
/AE iso (r) for different values of a, we gain insight in the re- 
lation between the ellipticity of the galaxies and their dark 
matter haloes. Note that as an alternative, we could weight 
Equation ([5]) with the lens ellipticity. 

It is useful to assess the signal-to-noise we expect to 
obtain for the shear anisotropy measurement compared to 
the signal-to-noisc of the tangential shear itself. For this 
purpose, we write Equation © in its most basic form: 

AS m odcl (r) = AE iso (r)[l + /cos(2A0)], (11) 

which has the following solution for the anisotropic part: 



/AEi, 



cos(2A6> t ) 
^ i m i cos 2 (2A0 1 ) 



(12) 



If the dark matter halo is described by a s i ngular isother- 
mal ellipsoid (SIE; see IMandelbaum et al.l l2006al) . and if 
the galaxy is perfectly aligned with the halo, we find 
/ = e h/2. Hence the anisotropic signal is a factor e/j/2 
lower than the isotropic signal. To assess the relative 
size of the error of /AEj SO compared to £; so , we insert 
Equation (fTTj) into Equation (|12[) , define a new weight w~l = 
Wi cos 2 (2A#i), and determine the error using c/A£ iso = 
VvEi^- Since wi and cos 2 (2A6>i) are uncorrelatcd, it 
follows that crj As . so = l/\/Y^i Wi(cos 2 (2A6)) = V2(TAS iso , 
with <TAs i30 = l/ySi w i t ne error on A£; so . Hence the er- 
ror of /AEi so is a factor \[2 larger than the error of A£j so . 
Consequently, the signal-to-noise of the anisotropic part of 
the lensing signal, (S/N) an i, is related to the signal-to-noisc 
of the isotropic part, (S/N)i so , as: 



< S / N >»' " (&) (S/N) ' 



(13) 



In the best-case scenario, the expected signal-to-noise of the 
shear anisotropy is an order of magnitude lower than the 
signal-to-noise of the azimuthally averaged shear. Applying 
the correction to remove systematic contributions increases 
the errors of the shear anisotropy by another factor of s/2. If 
the dark matter is described by an elliptical NFW, the sig- 
nal decreases rapidly with increasing separation (see Figure 
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2 of Mandclbaum et al. 2006a), and is only larger than 
the SIE signal on very small scales. If no redshift infor- 
mation is available for the lenses, the rapid decline of the 
shear anisotropy is particularly disadvantageous as the sig- 
nal can only be averaged as a function of angular separa- 
tion. Consequently, the anisotropy signal is smeared out, 
making it harder to detect. Finally, if the galaxy and the 
halo are misaligned, the signal decreases even further. These 
considerations show that we need very large lens samples 
to achieve sufficient signal-to-noise to enable a detection, 
and it motivates our choice to select broad lens samples. 

2.4. Contamination correction 

A fraction of our source galaxies are physically asso- 
ciated with the lenses. They cannot be removed from the 
source sample because we lack rcdshifts. Since these galax- 
ies are not lensed, but are included in calculating the av- 
erage lensing signal, they dilute the signal. To correct for 
this dilution, we boost the lensing signal with a boost fac- 
tor, i.e. the excess source galaxy density ratio around the 
lenses, 1 + f C g{r)- This is the ratio of the local total (satel- 
lites + source galaxies) number density and the average 
source galaxy number density. This correction assumes that 
the satellite galaxies are randomly oriented. If the satellites 
are preferentially radially aligned to the lens, the contam- 
ination correction for the azimuthally averaged tangential 
shear will be too low. If the radial alignment of the physi- 
cally associated galaxies has an azimuthal dependence, the 
shear anisotropy can either be biased high or low. 

This type of intrinsic alignment has been stud- 
ied with seemingly different results; some authors (e.g. 
lAgustsson fc Brainerdll2006t iFaltenbacher et al.ll2007l) who 
determined the galaxy orientation using the isophotal po- 
sition a ngles, have observed a stronger alignment than oth- 
ers (e.g. Iffirata et alj|2004t IMandelbaum et al.ll2005al) who 
used galaxy momen ts . T his discrep a ncy w as attributed by 
ISiverd et all ([2009D and lHao et all ([201 it ) to the different 
definitions of the position angle of a galaxy; the favoured 
explanation is that light from the central galaxy contami- 
nates the light from the satellites, which affects the isopho- 
tal position angle more than the galaxy moments one. As 
we measure the shapes of source galaxies using galaxy mo- 
ments, we expect that intrinsic alignment has a minor im- 
pact at most and can be ignored. 

To study whether the distribution of source galaxies has 
an azimuthal dependence, we perform the analysis sepa- 
rately using the galaxies residing within 45 degrees of the 
major axis, and within 45 degrees of the minor axis. On 
small scales, the extended light of bright lenses leads to er- 
roneous sky background estimates, which causes a local de- 
ficiency in the source number density. This deficiency is dif- 
ferent along the major axis and minor axis, which could bias 
the correction we make to account for physically associated 
galaxies in the source sample. To determine which scales are 
affected, we study the source number density around galax- 
ies as a function of their brightness and ellipticity. The re- 
sults are shown in Appcndix[Dj For galaxies with m r i < 19, 
we find a larger deficiency along the major axis on projected 
scales smaller than 0.2 arcmin; for galaxies with m r ' > 19, 
the deficiency is larger on projected scales smaller than 0.1 
arcmin. Therefore, we only use scales larger than 0.1 arcmin 
for galaxies with m r > > 19, and scales larger than 0.2 ar- 
cmin for galaxies with m r i < 19. The overdensities around 
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Fig. 4. Excess source galaxy density ratio as a function of 
projected distance to the lenses. The green squares (blue 
triangles) indicate the excess density ratio measured using 
sources within 45 degrees of the major (minor) axis. The 
arrows indicate the location of the virial radius at the mean 
redshift of the lenses. We find that the excess density ra- 
tio along the major axis is higher than along the minor 
axis, most noticeably for the 'red' sample. Please note the 
different scales of the vertical axes. 

the lens samples are shown in Figure [4] We find that the 
source sample is only mildly contaminated by physically 
associated galaxies, as the overdensities reach a maximum 
excess of only 30% for the 'red' lenses at the smallest pro- 
jected separations. The excess source galaxy density ratio 
is a few percent larger along the major axis than along the 
minor axis within the virial radii of the lens samples, most 
noticeably for the 'red' lens sample. 

The measured anisotropy is caused by two effects^: 
anisotropic magnification, and the presence of physically 
associated sources that are anisotropically distributed. As 
we lack redshifts for our galaxies, we cannot disentangle the 
two effects. However, we estimate the impact of anisotropic 
magnification for the lens samples in Appendix [Ej and find 
that even in the case where the galaxy and the dark matter 
halo are perfectly aligned, the effect is small. We conclude 
therefore that the observed anisotropy is the result of the 
anisotropy of the distribution of satellite galaxies. 

We correct the tangential shear in the major and minor 
axis quadrant for the contamination by satellites by mul- 
tiplying with their respective excess galaxy density ratio, 
before we measure the shear ratios. To calculate the cor- 
rection of (/ — As), we observe how /A£i so (r) changes in 
the presence of physically associated galaxies in the source 

4 Another effect is mentioned in IMan delbau m et all ([2006a) 
that could cause an anisotropic source density ratio: additional 
lensing by foreground galaxies. We estimate that this has a neg- 
ligible impact because the number of additional massive fore- 
ground galaxies is small due to our lens sample selection. 
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Fig. 5. Mean angle between the location of the satellites 
and the major axis of the lens galaxy as a function of pro- 
jected separation. The black triangles, purple diamonds and 
blue squares indicate the results for the 'all', 'red' and 'blue' 
lens sample. The arrows on the horizontal axis indicate the 
location of the virial radii at the mean redshift of the lenses, 
and correspond to 150 kpc, 280 kpc and 192 kpc for the 
'all', 'red' and 'blue' lens samples, respectively. The satel- 
lite galaxies preferentially reside near the major axis of the 
lenses. 

sample that are anisotropically distributed. Rather than 
Equation ([7]). the quantity we actually measure is 



/A£ iso (r) = - 



-V 



^AE.,e° t cos(2Afl. ; ) 
l + / cg (r,A0) ; 



iV = 2^m J e^cos 2 (2A0. i ), 



(14) 



with 1 + / cg (r, A#) the azimuthally varying excess galaxy 

density ratio, and ASj SO the unboosted lensing signal. We 
assume that 1 + / cg (V, A6*) has a similar azimuthal depen- 
dence as the shear, and can be described by 

1 + / cg (r, A0) = N iso (r) + 2N Ae (r)e« cos(2A0), (15) 

with a the exponent of the ellipticity used to weigh the 
shear measurement, iVj SO the azimuthally averaged boost 
factor and Nab the amplitude of the anisotropy. Using a 
Taylor expansion, we find that to first order 



f(r) = fir) + / eff (r), 



(16) 



with /eff(r) = AN A e(r)/N [so (r). To determine f e s{r), we 
measure both the angle-averaged boost factor, Ni so (r) = 
7Vls/-ZVlr, where ATls denotes the number of lens-source 
pairs and -/Vlr the number of pairs of lenses with ran- 
dom sources, and the azimuthally varying part, ^ A e(f) = 
Sls e g cos(2A6*)/-/Vlr. For the adopted model of the ex- 
cess galaxy density ratio this gives Ni so (r) = (1 + f cg (r))Ae, 



which is averaged over the angle, and £a0 = Ae(r)e 2 g ° 
These measurements are combined to give 



/eff(r) = A 



Ue(r) 



(l + /c g (r)>A*(e**>' 



(17) 



We determine the average value of f e g (r) within the virial 
radius, and add it to (/ — /4s). The values are tabulated 
in Table El Note that a similar correction is applied in 
iMandelbaum et al.l (|2006aD . 

To compare the anisotropy of the distribution of satel- 
lites to the literature, we now assume that at a narrow ra- 
dial range the excess galaxy density ratio can be described 



byl + /c: 



Nab cos(2A#). We fit this to the excess 



density ratio in the major and minor axis quadrants, sepa- 
rately for each radial bin. We use these fits to compute (8) , 
the mean angle between the location of the satellites and 
the major axis of the central galaxy, using 



(0) 



r 

Jo 



tt/2 



jZ /2 def CE {e) 



(18) 



In Figure [SJ we show (8) as a function of projected separa- 
tion for the three lens samples. 

Wc find that satellite galaxies preferentially reside near 
the major axis of the lenses, most strongly for the 'red' 
lenses. We determine the weighted mean of (8) within the 
virial radius, and find (8) = 43.7° ±0.3° for the 'all' sample, 
(0) = 41.7°±0.5° for the 'red' sample and (8) = 42.0°±1.4° 
for 'blue' sample. Additionally, for the 'red' lenses we find 
that (8) becomes more isotropic at larger projected separa- 
tions. It is useful to compare our results to previous stud- 
ies, that are based on simulations (e.g. Sales et all 20071 : 
iFaltenbacher et alj[200l lAgustsson fc Brainerdll201oh and 
observations fe.g. iBrainerdl 120051: lAgustsson k, Brainerdl 
I 2006L l201Ct IFaltenbacher et al.l 120071: iBailin et al.l 120081: 
iNierenberg et al.l 20111 ). In these works, (8) is found to be 
in the range between 41° and 43° for red central galaxies, 
whilst no anisotropy is observed for blue central galaxies. 
We can only make a useful comparison for the 'red' lens 
sample, as this sample is comparable to previously stud- 
ied red galaxy samples (i.e. predominantly containing red 
early-type galaxies, the majority of them expected to be 
centrals based on their luminosity distribution). We find 
that the constraints agree well. For the 'blue' and 'all' 
sample, we cannot make a comparison to previous work 
as these samples contain a mixture of early-type and late- 
type galaxies, and a fair fraction of them is expected to be 
a satellite of a larger system. The constraints we obtained 
are still interesting, however, as similar selection criteria 
can be applied to simulations, and the results compared. 

2.5. Virial masses and radii 

To determine to which projected separations the dark 
matter haloes of the galaxies dominate the lensing signal, 
we estimate the average halo size of each lens sample. For 
this purpose we model the azimuthally averaged tangential 
shear (after applying the contamination corrections) with 
an NFW profile, and fit for the mass. The NFW density 
profile is given by 



p(r) 



ScPc 



(r/r s )(l + r/r s y 



(19) 
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with 5 C the characteristic overdensity of the halo, p c the 
critical density for closure of the universe, and r s = 
r 20o/cNFW the scale radius, with Cnfw the concentration 
parameter. We ad opt the mass-concentration relation from 
IDuffv et all (120081) 



cnfw — J 



71 



M 



200 



2 X IQ^h^Mr. 



-0.084 



(l + z)-°- 47 , (20) 



Table 2. Shear ratios for the lens samples 



Sample 


U//mm) 


(i/C) 


All 
Red 
Blue 


i ic:+0.10 
- L - 1J -0.09 

u - y,:) -0.09 

1 16+ 019 


0.87 ±0.09 



which is based on numerical simulations using the best 
fit parameters of the WMAP5 cosmology. M 2 oo is defined 
as the mass inside a sphere with radius ^200, the radius 
inside of which the density is 200 times the critical density 
p c . We calculate the tangential shear profile using the 
analytical expressions pro vided by iBartelmannl (|l996h and 
IWright k Brainerdl (j2000ft . We fit the NFW profile between 
50 and 500 kpc at the mean lens redshift; closer to the lens 
the lensing signal might be contaminated by lens light, 
and at larger separations neighbouring structures bias the 
lensing signal high. The best fit Af2oo, r 200 and r s are given 
in Table [T] Note that in general, the best fit masses are 
lower than the mean halo mass because the shear of NFW 
profiles does not scale linearly with mass , and the distri- 
butio n of the halo masses i s not un i form (iTasitsiomi et al.l 
20041: iMandelbaum et all 120054 ICacciato et al.l l2009t 



Leauthaud et alj|2012t Ivan Uitert et all 1201 ll ). The result- 
ing uncertainty in the actual mass is not important here 
as we are mainly interested in the extent of the haloes, 
which is affected less (an increase of 30% in mass leads to 
an increase of only 10% in size). 



3. Shear ratio 

In this section we present the measurements of the 
ensemble-averaged ratio of the tangential shear along the 
major and minor axis of the lenses. This is a basic indica- 
tor of the presence of anisotropies in the lensing signal. Wc 
note that the shear ratio is not an optimal estimator as the 
weight is simply a step function, and does not depend on 
the elli pticity of the g alaxy . It enables, however, a compar- 
ison to iParker et al.l (|2007|) . Furthermore, we will use the 
shear ratio to examine how PSF residual systematics in the 
shape catalogues affect the anisotropy (Section l3~Tj) . 

For all elliptical non-power law profiles, the shear ra- 
tio varies as a function of distance to the lens. This radial 
dependence differs for diffe rent dark matter density profiles 
(|Mandelbaum et"a l. 2006a). Hence to obtain constraints on 
the halo ellipticity of the dark matter, we have to adopt a 
parti cular density profile . To compare our results to those 
from IParker et al.l (|2007| ) , we first assume that the density 
profile follows an SIE profile on small scales. In that case, 
the shear ratio is constant, and we determine the average 
and the 68% confidence limits as detailed in Appendix ICl 

In Figure O we show the average tangential shear along 
the major and minor axis, the average cross shear in the 
quadrants that are rotated by 45 degrees, and the inverse 
of the shear ratios / mm and /£% r . The tangential shear and 
the cross shear have been multiplied with the projected 
separation in arcmin, to enhance the visibility of the mea- 
surements on large scales where the signal is close to zero 
and the error bars are small. We sh ow the inverse of th e 
ratios following the definition used in IParker et al.l (|2007t ) . 
We do not observe a clear signature for an alignment or 



anti-alignment between the lenses and their dark matter 
haloes. Furthermore, we find that on small scales (<1 ar- 
cmin), / mm and f™™ are consistent, which suggests that 
the systematics present on these scales are smaller than 
the measurement errors. On larger scales, the difference is 
larger, which underlines the importance of applying the cor- 
rections to remove systematic contributions. The correction 
is largest for the 'all' lens sample, because its lensing sig- 
nal is smallest and therefore most susceptible to systematic 
contributions. We determine the average shear ratio within 
the virial radius at the mean lens redshift, and show the 
resu lts in Tabled 

IParker et all (|2007| ) used 22 square degrees of the 
CFHTLS to measure the shapes of ~ 2 x 10 5 lenses, se- 
lected with a brightness cut of 19 < i' < 22. Their lens 
sample consisted of a mixture of early-type and late-type 
galaxies with a median redshift of 0.4. The shear ratio was 
determined using measurements out to 70 arcsec (corre- 
sponding to 250 hr 1 kpc at z = 0.4), with a best-fit value 
of (l// mm ) = 0.76 ± 0.10. Excluding the round lenses with 
e < 0.15, the be st-fit ratio is (1//J = 0.56 ± 0.13. The 
lens sample from IParker et all ( 2007 ) can be best compared 
to our 'all' sample; comparing the relative number of early- 
/late-types in both sampl es using the CFHTLS W l photo- 
metric redshift catalogue ([Hildebrandt et aT]|2011l ) , we find 
they are similar. Also, the average mass of the lenses are 
comparable. Fitting the shear ratio on the same physical 
scale, we find (l//™"> = 0.98 ± 0.08 for the 'all' sample, 
which is ^2cr larger than IParker et ahl (|2007l ). Excluding 
lenses with e < 0.15, wc find (l//™m) = 0.95±g;^ s which 
is even almost 3cr apart. Since the lens samples are com- 
parable, this is most likely the resul t of differences in the 
analysis. Firstly, IParker et ahl (|2007l) do not apply a cor- 
rection for systematic contributions. However, systematic 
shear only tends to increase l// m m; if systematics were 
present, the discrepanc y would be even la rger. Secondly, 
it is not clear whether IParker et al.l (|2007l ) accounted for 
the non-gaussianity of the ratio of two gaussian distributed 
variables in determining the shear ratio; this is particularly 
important when the signal-to-noise of the lensing measure- 
ments is not very high. Generally, accounting for the non- 
gaussianity increases the positive error bar of the shear ra- 
tio, and decreases the negative one. This could bring their 
result closer to ours. Finally, it is not described how the av- 
erage ratio was determined. These differences could explain 
the discrepancy between the results. 



3.1. Imperfect PSF correction 

To measure the ellipticities of galaxies, we have to cor- 
rect their observed shapes for smearing by the PSF. The 
precision of the PSF correction is limited, which is mainl y 
due to the inaccuracy of the PSF model (lHoekstrall200l . 
Hence, residual PSF patterns may still be present in the 
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Fig. 6. Lensing signal multiplied with the projected separation in arcmin as a function of angular distance from the 
lens, for the 'all' lens sample (left-hand panels), the 'red' lens sample (middle panels) and the 'blue' lens sample (right- 
hand panels). In the top panels, the green squares (blue triangles) show the average rAE along the major (minor) axis 
(quadrants B (A) in Figure [3]). The dashed lines indicate the best fit NFW profile times the projected separation, fitted 
to the azimuthally averaged lensing signal on scales between 50 and 500 kpc using the mean lens redshift. In the middle 
panel, the green squares (blue triangles) show the cross shear signal averaged in quadrant D (C) of Figure[31 In the bottom 
panels, l// mm and l//^ r are shown by the red squares and black triangles, respectively. The dotted lines indicate the 
virial radius from the best-fit NFW profiles. The shear ratio does not provide clear signs for the alignment between 
galaxies and their dark matter haloes. 



shape catalogues. These residuals affect both the ellipticity 
estimates of the lens and the source galaxies, albeit with 
a different amount. Lens galaxies are typically large and 
bright, while source galaxies are small and faint, and hence 
harder to correct for. Regardless of that, PSF residuals 
tend to align the lens and source galaxies. If not accounted 
for, it could add a false anti-alignment signal to the shear 
anisotropy measurement (see Hoekstra et al. 2004) . 

We correct for PSF residual systematics in the cata- 
logues by subtracting the cross shear signal in the quad- 
rants that are rotated by 45 degrees with respect to the 
major and minor axes (7 K) c-d and / 45 A iso (r) in and 
(/ — /45)i respectively). To quantify how much PSF resid- 



uals actually contribute to these correction terms, and test 
whether they are properly removed, we introduce on pur- 
pose an additional bias in the PSF correction, and recalcu- 
late the shapes of the galaxies. Usually, the ellipticities of 
galaxies in the KSB method are computed as follows: 



f 



e- (1 + 6) x 



(21) 



with P 1 the shear polarisability, P sm the smear suscep tibil- 
ity tensor, and e the polarizations (jKaiser et alJll995l) . The 
starred quantities are determined using the PSF stars. The 
bias b is normally equal to zero, but to mimic an imper- 
fect PSF correction we set it to —0.05, and recalculate the 
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Fig. 7. Difference between the original and the PSF biased 
shear ratios (l//" r ) - (1 //,„"') bias (black triangles) and 
(l//mm) - (l//mm)bias (red squares) as a function of pro- 
jected separation from the lens for the three lens samples. 
We find that the PSF residuals are properly removed from 
the corrected ratio l//™ r , as the difference is consistent 
with zero on all scales. For the uncorrected ratio l// m m, 
the difference is negative and decreases with projected sep- 
aration. This result shows that the cross term effectively 
removes PSF residuals in the shear ratio estimators. 



shapes of all galaxies. We create new random shear cata- 
logues, and repeat the analysis using these biased shapes. 
We show the difference between the original and the biased 
shear ratios of the lens samples in Figure [7J 

We find that the difference of the shear ratios that 
are determined using the original and the PSF biased cat- 
alogues is consistent with zero on all scales for l//£^ r , the 
shear ratio estimator that is corrected with the cross shear 
terms. For the uncorrected shear ratio estimator, l// mm , 
we find that the difference is consistent with zero on small 
scales, but turns negative for projected separations larger 
than a few arcmins. This shows that if PSF residuals are 
still present in the shape catalogues, it affects l// mm , but 
not 1 / /m°m ■ Hence we conclude that PSF residuals are prop- 
erly accounted for using the cross shear signal. 



4. Impact of multiple lenses 

More than one lens may contribute to the shearing of 
a single source galaxy. Furthermore, some of the lenses are 
lensed themselves. In this section, we estimate the impact 
of these multiple lensing events on the halo ellipticity mea- 
surements. We also study the impact of the clustering of 
the lenses, and the correlation between their shapes, on the 
shear anisotropy. 



4.1. Multiple deflections 

Some foreground galaxies in our data lens both the 
lenses from the lens samples and the source galaxies. We 
denote these foreground galaxies with L2, and our selected 
lenses with LI. The impact of these 'multiple deflections' 
on the halo ellipticity me asurements were first discussed in 
iHowell k Brainerdl ((2010D . who found that it adds a strong 
false anti-alignment signal to the shear anisotropy measure- 
ments. Multiple deflections affect the halo shape measure- 
ment in three ways. Firstly, the orientation of the lens light 
changes, leading to a misalignment if the light was aligned 
with the halo. Consequently, the lensing signal is averaged 
in quadrants that are rotated with respect to the unlensed 
ones, causing a reduction of the shear anisotropy. This ef- 
fect is only important if the projected separation between 
L2 and LI is small, as only those configurations lead to 
significant changes in the ellipticity of LI. Secondly, the 
source galaxies experience shear not only from LI, but also 
from L2. Especially source galaxies close to L2 are affected. 
Finally, the lensing of L2 changes the observed positions of 
LI and the sources. We ignore this third effect as the im- 
pact is negligible. 

In the presence of L2, Equation (Ql changes to 



Jmm{r) 37r/4+(5e ,7tt/4+<50 - M ' 

L/4+se ItdS + h./i+se ^ dd 



(22) 



where the integration is performed over quadrants that are 
rotated by 58, the change of the position angle of LI caused 
by the lensing of L2. j t is the sum of the shear of LI and 
L2 at the location of a source galaxy. The 71-component 
is given by 71 = 74,1,1 cos(20) + j t ,L2 cos(20), with j t , L i 
and 7t,L2 the tangential shear of LI and L2, 9 the angle 
between the source galaxy and LI and <f> the angle between 
the source galaxy and L2. Hence the signal that is measured 
is given by 



li 



;COs(2^)] x cos(26»)- 



[j t ,Ll C0S(26») + J t ,L2< 

foil sin(2<9) + 7t,i 2 sin(2^)] x sin(20). 



(23) 



The change of the equations for /™m /A£i so (0 and 
/45ASj SO (r) in the presence of L2 can be derived in a simi- 
lar way. 

To obtain an intuitive understanding of the impact of 
multiple deflections on the halo ellipticity measurements, 
we compute the change of the shear anisotropy of a single 
lens in the presence of an additional foreground galaxy us- 
ing simple idealised simulations. These simulations, which 
are discussed in Appendix [I] suggest that multiple deflec- 
tions mainly affect the shear anisotropy of round (e < 0.15) 
lens galaxies, and at large projected separations. To confirm 
these findings, we create a large set of simulated image cat- 
alogues to obtain a rough estimate of the impact. 

For the simulated catalogues we adopt an image size 
of 30x30 arcmins. We randomly assign positions to 10 000 
galaxies (approximately the galaxy number density of the 
RCS2). Redshifts are assigned to each background galaxy 
by drawing from the redshift distribution of the RCS2 
source galaxies. The background galaxies are intrinsically 
round, and not convolved with a PSF, to avoid introduc- 
ing unnecessary sources of noise. We insert 50 lenses at a 
typical lens redshift of z = 0.4, to which we assign an el- 
lipticity e g with a random value between and 0.4 and a 
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random position angle. The ellipticity of the dark matter 
halo, eu, is proportional to the ellipticity of the galaxy via 
e h = fh x G g . We use a fixed value for /h = 1.0, but we also 
test the impact of loosening this assumption. Each lens is 
modeled with an SIE prof ile, and the induced shear o n each 
source is computed with (|Mandelbaum et aT1l2006a[) 



li 



Ana 2 DiDis 1 
c 2 D, 2r 



1 + y cos(2(9) 



(24) 



with eh the ellipticity of the dark matter halo, and a the 
velocity dispersion of the lens which we have set to 200 km 
s _1 . We create 20 sets of 500 catalogues, and determine 
the mean lcnsing signal and the scatter between the 
simulations sets. This enables us to assess the significance 
of potential trends. 

To study the impact of multiple deflections, we would 
ideally assign velocity dispersions to all galaxies that reside 
in front of the lenses following their velocity dispersion 
distribution, and use them to compute the shear on the 
lenses and sources. The galaxies that reside behind the 
lenses only introduce noise, and can be ignored. This is 
computationally expensive as there are many galaxies at 
lower redshifts. However, the majority of the foreground 
galaxies are not massive, and are close in redshift to the 
lenses (resulting in small lensing efficiencies), so their 
contribution to multiple deflections is negligible. 

For computational speed-up, we therefore define a 
smaller number of foreground galaxies. We choose their 
velocity dispersions and redshifts such that the impact 
of multiple deflections is comparable to what is expected 
using all the foreground galaxietQ. On these grounds, we 
randomly insert a second set of 500 round lenses with a 
truncated isothermal sphere (TIS) profile with a velocity 
dispersion of 100 km s _1 and truncation radius of 150 
arcsecs, located at a redshift of 0.1. Both the source 
galaxies and the LI lenses are lensed by the L2 lenses, and 
we change their ellipticities accordingly. Then we measure 
the shear anisotropy around the LI lenses as we would 
in observations, i.e. using the 'observed' ellipticities. We 
show the anisotropy of the lensing signal in panel (a) of 
Figure |5J confirming the predicted trends from Appendix 
[Fl multiple deflections lead to a reduction of the shear 
anisotropy, with a magnitude that increases for larger 
separations to the lens. On small scales, the reduction of 
the shear anisotropy is larger for the shear estimators that 
have not been corrected for systematic contributions (/ 
and / mm ) than for the corrected ones ((/ — As) and /m° m r ) 
(see Figure |R2p). As long as the separation between the 
sources and LI is small, the additional shear from L2 is 
relatively constant and hence efficiently removed using the 
cross terms. The correction does not work on larger scales 
as the additional shear from L2 varies spatially. Note that 



5 Using the photometric redshift catalogu e of the CFHTLS 
"Deep Survey" fields from lllbert et al] (120061 ). we find that ap- 
proximately 30% of all galaxies in the RCS2 have a redshift 
z < 0.4, with a mean of z = 0.35. The lensing efficiency at 
z = 0.1 is about 6 times larger than at z — 0.35 for a source 
at z = 0.4, hence we can reduce the number of L2 lenses by a 
factor of six when we place them at z — 0.1. Note that we ignore 
the increase of the lensing efficiency by a factor of two for the 
source galaxies when the L2 lenses are placed at z = 0.1, which 
increases the impact of multiple deflections. Additionally, the 
use of a smaller number of more efficient lenses also leads to an 
increase of the impact of multiple deflections. 



the reduction of the shear anisotropy is smaller for (/ — /4s) 
than for /m™ 1 , because the signal of the former is weighted 
with the ellipticity of the lenses; the most elliptical lenses 
are less affected by multiple deflections (see Appendix IF)) . 

Based on Figure IF. 2b in Appendix we expect that 
the impact of multiple deflections is reduced if we exclude 
the roundest lenses. Therefore, we repeat the simulations, 
excluding lenses with an observed (rather than intrinsic) 
ellipticity e < 0.05 as is done in the measurements on the 
real data. We show the results in panel (b) of Figure [51 
Excluding the roundest lenses significantly reduces the 
impact of multiple deflections for / mm and /m°m- The 
improvement for / and (/ — /4s) is minor, as the roundest 
lenses are already downweighted in this measurement. 
There is some residual signal left on the largest scales, but 
to constrain halo shapes in real data we only use measure- 
ments on small scales. Therefore, we find it unlikely that 
multiple deflections strongly biases the lensing anisotropy. 

These conclusions do, however, depend on our as- 
sumptions. The number and the masses of the L2 lenses 
is most critical. In our simulations, we have assumed 
an average velocity dispersion of 100 km s" 1 for the 
L2 lenses. However, more massive L2 lenses contribute 
more to multiple deflections as their shear patterns affect 
larger patches of the sky. Most of the massive L2 lenses 
reside relatively close in redshift to the LI lenses, and 
consequently their lensing efficiency is small. Therefore, we 
do not expect that assuming a constant velocity dispersion 
for the L2 lenses rather than drawing them from the 
distribution has a large impact on the results. 

A second important assumption is that we have 
modeled the density distribution of the LI galaxy with 
an SIS. This is a reasonably accurate description of a 
galaxy density profile at small scales, but not on large 
scales. For NFW profiles, which are more appropriate, the 
lensing anisotropy declines strongly as a function of radius 
(jMandelbaum et al"1l2006al) . A small reduction by multiple 
deflections of an already small signal has a relatively 
larger impact. One way to account for this is to give larger 
weights to the measurements on small scales. Alternatively, 
the shear anisotropy could be fit using measurements on 
small scales only. 

We also assumed that the relation between the ellip- 
ticity of the galaxy and the ellipticity of the dark matter 
halo is linear. To study the impact this might have on the 
anisotropy of the lensing signal, we loosen the assumption. 
We assign a random value to e/j between and 0.4, but 
we keep the position angles of the dark matter and the 
galaxy perfectly aligned. Not surprisingly, we find that 
the shear ratios / mm and f™™ arc unchanged, as the 
average halo ellipticity does not change. Without multiple 
deflections, we find that (/ — / 45 ) is reduced by 25% to 
0.1875. In this case, the anisotropy of the lensing signal 
is no longer proportional to e 2 , but to e/j x e g : the factor 
e/i comes from Equation (|24l) . and the factor e g from the 
weight in Equation ([7]). When we average over a flat lens 
ellipticity distribution, the lensing signal in the original 
case is oc J de g e 2 , which becomes cx (e^) J de g e g if e/,. 
is randomly assigned. The second integral is 25% smaller 
than the original. Hence the decrease of (/ — /4s) results 
from no longer giving a larger weight to the galaxies with 
large dark matter halo ellipticities in the measurement. 
If we include multiple deflections, the relative decrease of 
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Fig. 8. (a) Anisotropy of the lensing signal in the simulations in the presence of additional foreground galaxies that lens 
both the lenses and the sources. In the top panel, we show (/ — /4s) and / with the filled black and open green diamonds, 
computed with the exponent in the weight a = 1.0. In the lower panel, we show f^™ (/mm) with the filled black diamonds 
(open green diamonds). The dashed lines show the signal in the absence of multiple deflections. The uncorrected shear 
estimators / and / m m are increasingly reduced for larger distances to the lens. The reduction of the systematic shear 
corrected estimators is clearly smaller on scales < 100 arcsec, but the difference decreases at larger separations as the 
systematic shear is not constant anymore and the correction therefore inaccurate, (b) Anisotropy of the lensing signal in 
the simulations in the presence of additional foreground galaxies, excluding lens galaxies with e < 0.05 from the analysis. 
The exclusion of round lenses is found to reduce the impact of multiple deflections, particularly for / mm and f™m as 
these measurements are not weighed with the lens ellipticity. 



(/ — /4s) is similar as in Figure [8^. 

Our simulations show that on small scales (<1 arcmin), 
the impact of multiple deflections is a few percent at most. 
This result is robust to changes in the simulation set-up, 
which shows that it is very unlikely that multiple deflec- 
tions affect the shear anisotropy measurements on these 
scales. On larger scales, the impact of multiple deflection 
is more uncertain, and could have an important effect 
on weak lensing studies that investigate the alignment 
between galaxies and the large-scale structure. Ultimately, 
realistic numerical simulations should be used to accurately 
assess and model the impact multiple deflections have 
as a function of the selection of the lens sample. Such 
simulations are also essential to investigate the impact 
of large-scale structure and voids, and test whether their 
impact on the shear anistropy on scales close to the lens 
is small as expected. This is necessary to improve the 
precision of the gravitational lensing constraints on the 
shapes of dark matter haloes with data from upcoming 
sur veys. 

iHowell fc Brainerdl f|2010h report a significant and 
strong decrease in / mm with projected lens-source 
separation. Howe v er, f or computational speed-up, 
IHowell fc Brainerdl (|2010h only use those foreground 
galaxies that reside within 100 arcsec from the lens to 
calculate the multiple deflection signal. Figure fF2k shows 
that this leads to an overestimation of the reduction of the 
shear ratios. Furthermore, redshifts are assigned to galaxies 



whose positions have been observed in real data. This 
causes an unphysical correlation between the positions 
of foreground and background galaxies, and amplifies 
the chance of having a low-redshift foreground galaxy 
close to a lens, which could result in an overestimate of 
the effect. Finally, no correction for systematic shear is 
implemented, which would have reduced the impact as well. 



4.2. Clustering of lenses 

Many galaxies in the Universe reside in groups and clus- 
ters. Hence the shear we measure around a lens galaxy is 
the sum of the shear from the lens and neighbouring galax- 
ies. At small projected separations from the lens, the signal 
is dominated by the lens galaxy, but at larger separations 
the shear from the neighbours becomes increasingly impor- 
tant. This may affect the shear anisotropy measurements. 
To study this, we assign the positions of the lenses from 
the real lens samples to the simulated lenses. The lens el- 
lipticitics arc randomly drawn between and 0.3, and are 
modeled with an SIE profile. The positions of the galaxies 
with 22 < m r i < 24 in the RCS2 are assigned to the sim- 
ulated sources. The lenses are put at a redshift z = 0.4, 
and the source redshifts are randomly drawn from their 
redshift distribution. We assume that /h = 1, and that 
no foreground lenses are present as that would potentially 
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Fig. 9. (a) Anisotropy of the lensing signal in the simulations using the observed positions of the 'red' lens sample in 
the RCS2. The filled black diamonds (open green diamonds) show the shear anisotropy estimators (not) corrected for 
systematic shear. The cllipticities of the lenses have been randomly assigned. The lensing anisotropy declines slightly at 
scales > 20 arcsec, which is expected as clustering of lenses with random position angles leads to an isotropic signal at 
large scales. In panel (6), we add an additional ellipticity of 0.03 to the e\ component of the lenses to mimic intrinsic 
alignments. This increases the shear anisotropy, but also induces cross shear, so that the systematic shear corrected shear 
anisotropy estimators are unchanged. 



mix different effects. To obtain the errors, we determine 
the scatter between 20 random realizations. We show the 
results in Figure for the 'red' lens sample. 

We find that clustering of the red lenses slightly re- 
duces the lensing anisotropy on large scales. This can be 
easily understood: if galaxies with random ellipticities clus- 
ter, the resulting lensing signal around the lens ensemble 
at large scales becomes more isotropic. The impact of clus- 
tering on the shear anisotropy around the 'blue' and 'all' 
lenses is smaller, and is therefore not shown. We note that, 
additionally, dark matter haloes may be stripped in high- 
density environments, which could reduce the anisotropy 
signal as well. 

Numerical simulations suggest that dark matter haloes 
are a l igned, and point toward each other fe.g. | Splinter et al 
1997tlCroft k Met zlcr 20001 iHeavens et al.ll2000t iLee et al 



20081) . If the galaxies that reside in these haloes are prefer- 
entially aligned with the halo - a prerequisite for measuring 
halo shapes - the observed galaxy ellipticities are expected 
to show this alignment as well. In particular, the elliptici- 
ties of luminous red galaxies are believed to b e increasingly 
correlated with decre asing separations (e.g . Hirata et al.l 
I2007t lOkumura et"all 120091: Uoachimi et all 1201 ID . which 
can be understood in a framework in wh ich galaxies form 
in a linear tidal field (iBlazek et al.l [20TlT ). Similar studies 
for spiral galaxies show t he effect is considerably weaker 
(jMandelbaum et "all 1201 1( ). Intrinsic alignments also affect 
the shear anisotropy signal. 

To obtain a conservative estimate of the impact of in- 
trinsic alignments, we add an ellipticity of 0.03 to all the e\ 
components of the lenses. The result for the 'red' lenses is 



shown in Figure|9jD. We find that correlated lens ellipticities 
increases the shear anisotropy on large scales; the shear pat- 
tern of neighbouring lenses amplifies each other. However, 
correlated lens ellipticities also induces cross shear of an 
equal magnitude, hence the corrected shear anisotropy es- 
timator is unchanged. In reality, the intrinsic alignments 
are scale dependent, and the correlation decreases for larger 
separations. Furthermore, if the lens galaxies are pointing 
towards each other, the contributions of the shear along the 
major axis systematically add up, which could amplify the 
effect. This is difficult to model, however, but we will assess 
it with numerical simulations in a future work. 

The main conclusion from our simulations is that the 
impact of multiple deflections, clustering of lenses and the 
correlations of their cllipticities is small, and can be safely 
ignored for the lens selection in this work (i.e. particularly 
for the 'blue' and 'red' lens sample, which were selected to 
contain massive, elliptical, and low-redshift galaxies). More 
realistic studies using numerical simulations are required 
to quantify the effects more precisely. This is a crucial step 
in correctly interpreting the lensing anisotropy in future 
lensing data that are of higher precision. 

5. Halo ellipticity 

The simulations from the previous section indicate that 
multiple deflections and clustering have a small impact on 
the shear anisotropy, in particular for the 'blue' and 'red' 
lens samples where we have selected massive, elliptical, low- 
redshift lenses. Now we proceed with the actual measure- 
ments. We show (/ — /45)AS; so as a function of projected 
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Fig. 10. Anisotropic shear signal around the 'all' lenses (top panels), the 'red' lenses (middle panels) and the 'blue' 
lenses (bottom panels) . The signal has been multiplied with the projected separation in arcmin to enhance the visibility 
of the signal on large scales. The green triangles denote r/A£i so , the red squares r/4 5 ASi SO , and the black diamonds 
r(f — /45)A£j so . The vertical dotted lines indicate the virial radius at the mean lens redshift, the dashed lines show the 
best fit elliptical NFW profiles. Indicated on top of each column is the weight applied to the measurement. The average 
shear anisotropy of the 'all' and 'red' lenses is weakly positive within the virial radius, supporting the presence of triaxial 
dark matter haloes that are reasonably aligned with the galaxies. For the 'blue' lenses the signal is marginally negative. 
Note that the measurements shown here are not corrected for the anisotropic distribution of satellite galaxies. 



distance to the lens in Figure [TU1 We multiply the lensing 
signal with the projected separation in arcmin to enhance 
the visibility of the signal on large scales. Each row shows 
the measurements for one of the lens samples, whilst the ex- 
ponent of the galaxy ellipticities, a in Equation ([9]), differs 
between the columns. To quantify the shear anisotropy, we 
first determine the average value of (/ — f^) as detailed in 
Appendix [C] The best-fit values are summarized in Table 

m 

For the 'all' sample of lenses, we find that (/ — ^45) is 
marginally positive, independent of the applied weight. For 
the 'red' lenses, we find that (/ — f^) is consistent with 
zero if we take w oc e . For larger values of a, however, 
(/ — 745) turns marginally positive due to the two data 
points near the virial radius. Closer to the lens, the sig- 
nal is slightly negative. For the blue galaxies, the signal is 
marginally negative. 

Next, we fit an SIE and an elliptical NFW profile. For 
the SIE fit, we first determine the Einstein radius Te by 
fitting a singular isothermal sphere (SIS) profile to the az- 



imuthally averaged tangential shear measurements within 
the virial radius at the mean lens r edshift. To determine 
/ h , we fit (jMandelbaum et al.ll2006al) 

(/-/ 45 )AS(r) = ^pA; (25) 
8r 

the errors of /h are determined from the x 2_va l ues °f the 
fit. Since (,f — f^)/ fh = 0.2 5 on all scales for an SIE profile 
dMandelbaum et aLll2006aT) . we expect to find best-fit val- 
ues for /h that are four times larger than (7-/45) . The best- 
fit values are shown in Tabled We find that /h is consistent 
with, but not exactly, four times (/ — f^). The difference 
may be due to differences between the fitting methods; if 
we first fit an SIS profile, we correlate the azimuthally aver- 
aged tangential shear measurements, whilst if we determine 
(/ — /is) directly from the data, each data point is treated 
separately. The general trends, however, are consistent and 
the conclusions do not depend on how we fit the data. 

Currently, no analytical expression exists for the shear 
anisotropy of an elliptical NFW profile. Therefore, we use 
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Table 3. The best-fit values for the anisotropy of the galaxy-mass cross-correlations function, (/ — /4s), and the ratio of 
the dark matter halo ellipticity and the galaxy ellipticity, /h, for an SIE and an elliptical NFW profile. 

Sample a (jW) (/ - / 45 > ApE) A (NFW) 



All 


0.0 


1.3 ±0.6 x 10" 


-3 


0.19 ±0.10 


0.47 ±0.37 


All 


()..-) 


1.1 ±0.7 x 10" 


-3 


21 +0 11 
u - zi -o.io 


0.57 ±0.40 


All 


1.0 


0.8 ±0.8 x 10" 


-3 


0.23 ±0.12 


0.70 ± 0.46 


All 


1.5 


0.6 ± 1.0 x 10" 


-3 


0.26 ±0.15 


0.83 ±0.55 


All 


2.0 


0.4 ± 1.2 x lO- 


-3 


0.29 ±0.17 


0.97 ±0.65 


Red 


0.0 


ll. 9 ± 1.8 x 10 


-3 


0.13 ±0.15 


0.00 ±0.58 


Red 


0.5 


11.3 ± 2.1 x 10 


-3 


0.19 ±0.16 


0.05 ± 0.60 


Red 


1.0 


9.3 ±2.5 x 10" 


-3 


0.28 ±0.18 


0.25 ±0.70 


Red 


1.5 


7.2 ±3.1 x 10" 


-3 


0.40 ±0.22 


0.61 ±0.86 


Red 


2.0 


5.2 ±4.0 x 10" 


-3 


0.54 ±0.27 


1.09 ± 1.07 



0.96Z 
1.19j 



H0. S3 
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1.801 

2.121 
-0.19; 



-1.09 
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-0 14+ 1 12 

0.20±i^ 

0.87ii;g 
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1.82; 



Blue 


0.0 


1.5 ± 1.4 x 10" 
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-0 16 +0 18 


-0.56 ±0.68 


Blue 
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-0.25 ±0.19 


-0.75 ±0.70 


Blue 
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-3 


U.OO_g 22 


-1.01 ±0.81 


Blue 
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2.5 ±2.3 x 10" 


-3 


-0.45 ±0.26 


-1.24 ±0.96 


Blue 


2.0 


2.5 ±2.7 x 10" 


-3 


u - dJ -0.32 


-1.44 ± 1.17 



1.62 
1.65 
1.69 
1.72 
1.97 
2.03 
o f?7+2. 36 
'• Dl -2.44 

-3.06±2;95 



-1.247 
-1.62j 
-2.I7] 



the numerically integr ated values of //./>, and f^ /fh as a 
function of r/r s from iMandelbaum et al.l ((20063) (shown 
in Figure 2 of that paper), which have been kindly pro- 
vided by Rachel Mandclbaum. Since our lens galaxies span 
a broad range in redshifts, we first determine the redshift- 
averaged lensing model by integrating the elliptical NFW 
profiles over the redshift distribution of each lens sample 
(shown in Figure [IJ , and weigh each lens redshift bin with 
the lensing efficiency (Di s /D s ) that is averaged over the 
source redshift distribution. Note that this is an important 
correction; for the 'all' sample, the integrated profile results 
in about 50% larger values for (/ — /4s) compared to the 
profile computed using the mean lens redshift. For the 'red' 
and 'blue' sample, the difference is smaller because their 
redshift distributions are narrower. Also note that the SIE 
profiles do not have to be corrected, since the azimuthally 
averaged lensing signal and the anisotropic part are simi- 
larly affected when we integrate them over the lens redshift 
distribution because both scale as r _1 . The best-fit values 
of fh are therefore unaffected. 

The best fit values of fh for the elliptical NFW pro- 
files are less significant than (/ — /4s) for the same lensing 
measurements. The reason is that the elliptical NFW fit is 
very sensitive to the signal close to the lens, but not to the 
signal at larger separations. We find that for the 'all' and 
'red' sample, (/ — /4s) actually turns slightly negative close 
to the lens, rather than increasing strongly as would have 
been expected for an elliptical NFW profile that is aligned 
with the lens. Although this might be just caused by noise, 
it could also indicate that a single elliptical NFW profile 
does not describe the shear anisotropy signal well. 

Finally, we note that for the 'all' sample (/ — ^45) AE; so 
turns negative at projected separations >5 arcmin. A sim- 
ilar trend can be observed in Figure [6j where the inverse of 
the corrected shear ratio of the 'all' sample is slightly larger 
than unity. We cannot directly interpret this as the result of 
an ant i- alignment of galaxies with the large-scale structure, 
as we found in the previous section that multiple deflec- 



tions and the clustering of galaxies produce a similar trend 
at these scales, and we cannot disentangle the effects. To 
constrain the average halo cllipticities of galaxies we only 
use the lensing signal on scales <1 arcmin, however, where 
the effect of multiple deflections and clustering of galax- 
ies can be safely ignored, and a non-zero signal reflects an 
anisotropy of the projected gravitational potential. 

5.1. Environmental dependence 

To study whether the lensing anisotropy depends 
on the lens environment we measure the signal for the 
isolated and clustered lens sample. In Figure [11] we show 
(/ — /45)ASi SO for the w oc e 1,0 bin, which is the same 
weight as used in IMandelbaum et aLl (|2006al ) and hence 
enables a direct comparison. We determine (/ — f^) and 
fit /h for the elliptical density profiles, and show the results 
in Table [U 

The lensing anisotropy for the isolated 'all' lenses is 
positive, and the values of (/ — fa) ASi SO and fh are larger 
than those of the clustered sample by almost ~2er. For the 
'red' and the 'blue' lenses, we find that on small scales, 
the lensing anisotropy is more negative for the clustered 
sample. When we average the signals within the virial 
radius, or fit the elliptical density profiles, we find that 
this difference is not statistically significant. 



5.2. Interpretation 

The shear anisotropy measurements provide weak 
support that the average galaxy is preferentially aligned 
with its triaxial dark matter host. The significance of the 
detection for the 'all' sample does not depend on how we 
weigh the measurement with the observed galaxy elliptic- 
ity, which indicates that more elliptical galaxies do not 
reside in, or are better aligned with, more elliptical dark 
matter haloes. We find that the errors on (/ — f^) increase 
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Table 4. Similar to Table [3] for the isolated and clustered lenses. 



Sample 


a 


Environment 


(/ - fa) 


/h(SIE) 


/h(NFW) 


All 
All 


1.0 
1.0 


isolated 
clustered 


0.5ltoi 6 5 
0.18 ±0.13 


2.35 ± 1.00 
0.40 ±0.51 


0.9011;" 


Red 
Red 


1.0 
1.0 


isolated 
clustered 


0.11 ±0.25 
0.47 ±0.26 


0.49 ±0.97 
0.46 ± 1.05 


40+ 1 96 

1 QQ+1.90 

1 - zo -1.83 


Blue 
Blue 


1.0 
1.0 


isolated 
clustered 


-0.18 ±0.30 

-0 4S+ 34 


0.05 ± 1.06 
-1.54 ± 1.33 


0.38tJ;g 
-3.46± 3 3 ; 2 J 
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Clustered 
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Fig. 11. Anisotropic shear signal multiplied with the pro- 
jected separation in arcmin around the 'all' sample (top 
panels), the 'red' sample (middle panels) and the 'blue' 
sample (bottom panels), for the isolated lenses on the 
left-hand side, and the clustered lenses on the right-hand 
side (using w oc e 10 ). The green triangles denote /ASj so , 
the red squares /4sAEi so , and the black diamonds (/ — 
/45)AEj S0 . The vertical dotted lines indicate the virial ra- 
dius at the mean lens redshift, the dashed lines indicate the 
best fit elliptical NFW profiles. Please note the different 
scalings of the vertical axes. For the 'all' sample, the shear 
anisotropy is larger for the isolated sample; for the 'red' 
and 'blue' sample, the shear anisotropy is more negative on 
small scales for the clustered samples, but the differences 
are not significant when averaged within the virial radius. 



for larger values of the exponents of the weights, because 
the effective number of lenses decreases: round lenses 
barely contribute to the signal for the largest exponents. 
Furthermore, we find that the shear anisotropy signal of 
isolated galaxies is stronger than that of clustered galaxies: 
the difference is almost 2a. One possible explanation 
is that the dark matter haloes of isolated galaxies arc 
less subject to stripping, and may preserve their original 



shapes. Clustered galaxies, on the other hand, may lose 
more of their dark matter, particularly in the outer regions. 
An alternative explanation is that the fraction of lenses 
in the clustered sample that are satellites is larger. Since 
the host halo in which satellite galaxies are embedded 
dominates the lensing signal, and since the orientation of 
the major axis of the satellites and the major axis of the 
central galaxy are expected to not be strongly aligned, 
this would also lead to a reduction of the anisotropy of the 
lensing signal. 

The average shear anisotropy of both the 'red' and 
'blue' lenses is consistent with zero. We find no significant 
differences between the clustered and isolated sample, be- 
cause we lack precision. For the 'red' lenses the increase of 
the detection significance of (/ — f±$) to ^2a by increasing 
the exponent of the lens cllipticity in the weights could 
be interpreted as an indication that red galaxies with 
larger ellipticities preferentially reside in more elliptical 
dark matter haloes, and/or are better aligned with them. 
However, Figure [10] reveals that the detection is caused 
by the two measurement points near the virial radius, 
and not by the measurements at the smallest scales where 
the signal should be highest if the galaxies arc aligned 
with their triaxial dark matter haloes. This is reflected 
by the best-fit values of fh for the elliptical NFW profile, 
which are consistent with zero for all exponents. It seems 
therefore unlikely that the detection is caused by triaxial 
dark matter haloes. For the 'blue' lenses, we find that 
for the largest exponents the anisotropy is marginally 
negative, which is suggestive of an anti-alignment between 
the galaxy and the dark matter halo. 

The negative shear anisotropy at small projected sepa- 
rations for the 'red' lens sample occurs at the same scale 
where the contamination is highest. Since we observed in 
Section l2Tl that the distribution of satellites is anisotropic, 
the observed signal could also be caused by a radial 
alignment of physically associated galaxies in the source 
sample. We estimate the value of the average tangential 
intrinsic alignment that would produce the signal we 
observe in Appendix [G] We find that we require a value 
that is roughly te n times larger then the results from 
iHirata et al.l (|2004D . Conversely, the expected impact of 
intrinsic alignments is about a factor ten smaller than the 
signal we observe, hence it is unlikely that our measure- 
ments are significantly affected. Hence the negative shear 
anisotropy on small scales is unlikely caused by intrinsic 
alignments. 

Currently, weak lensing constraints on the ellipticity of 



17 



Edo van Uitert et al.: Constraints on the shapes of dark matter haloes 



dark m a tter haloes have bee n present ed in IHoekstra et al 
(2001, iMandelbaum et all (|2006al) and iParker et al 
2007 ). We already co mpared our results with those 
from IParker et ah! (|2007l) in Section [31 IMandelbaum et ahl 



(|2006air measured the weak lensing anisotropy around 
2 x 10 6 lenses with photometric redshifts in the SDSS. 
On average, they find //, = —0.06 ± 0.19 for the red 
galaxies, and fh = —1.1 ±0.6 for the blue galaxies on scales 
20-300 ft, _1 kpc for an SIE profile. Fitting an elliptical 
NFW profile yields f h = 0.60 ± 0.38 and f h = -lAt\ 7 Q 
for the red and blue lenses, respectively. Separating the 
lenses in luminosity bins, they find that the lensing signal 
is consistent with zero for most of their bins. Only for 
the brightest red lenses a detection is reported, also 
at the ^2 sigma level. A detailed comparison between 
the results is complicated due to differences in the lens 
selections. Our 'red' lens samples can be best compared 
with th e L5 and L6 red sample from IMandelbaum et ahl 
(l2006af). as these are most similar in absolute magnitudes; 
IMandelbaum et all ((2006al) find (/ - / 45 ) = 0.08 ± 0.08 
and (/ — f/t§) — 0.29 ± 0.12 for L5 and L6, which agrees 
well with the (/ - / 45 ) = 0.28 ± 0.18 that we found. For the 
elliptical NFW profile, IMandelbaum et al.l (|2006al) found 
f h = 0.4 ± 0.57 and f h = 1.7 ± 0.7 for the red L5 and L6 
sample, and the average of those values is roughly within 
lcr of our best fit value 0.20l}'g^. Our 'blue' lens sample 
covers a broad range in luminosity, and roughly corre- 
sponds to the blue lens b ins L3 to L5 ofpVlandelbaum et ahl 
(|2006al) . For these bins, IMandelbaum et ahl (|2006aD finds 
(/ - As) = -0.29i°;^, -0.36t°; 2 | and -0.27 ± 0.28 

respectively, which agrees well with —0.35^0 22- I n con ~ 
elusion, we find t h at our results and the results from 
IMandelbaum et ahl (|2006a|)_ are consistent. 

In an earlier work, iHoekstra et al.l (l2004h used 45.5 
deg 2 of the RCS (|Gladders fe Yed I2005T) to measure the 
lensing anisotropy around 1.2 x 10 5 lenses selected with 
19.5 < R c < 21. A -2cr detection of f h = 0.77^2? was 
obtained by fitting a TIS using a maximum likelihood 
method. This result appears to be signi fi cantly different 
from ours, and from IMandelbaum et ahl (|2006al) . In the 
latter, various reasons are given why the results could 
differ, which also apply to us: most importantly, the lens 
samples are very different, and the method of analysis 
differs. Possibly, the maximum likelihood method is a 
better estimator of the shear anisotropy, as it takes into 
account the positions and relative orientations of the 
lens galaxies. The result from a maximum likelihood fit, 
however, is difficult to interpret as it is not clear how 
multiple deflections affect the measurement, and because 
the radial dependence of the signal cannot be visualized, 
which makes it difficult to identify residual systematics. 



Converting the average shear anisotropy into the 
average projected ellipticity of dark matter haloes is 
complicated. The simulations in this work have shown that 
multiple deflections and clustering can have some impact 
on the lensing anisotropy. If we select the lens sample 
carefully as we did (in particular the 'blue' and 'red' lens 
sample that consist of isolated, elliptical, low-redshift, 
massive lenses), the impact of these complications is small. 
However, misalignments between the position angles of the 
galaxies and their dark matter haloes cause a reduction 
of the lensing anisotropy. Therefore, the observational 



constraints on the shear anisotropy from previous weak 
lensing analyses, as well as from this work, only provide 
lower limits on the average halo ellipticity. To constrain 
the halo ellipticity from the lensing anisotropy, we need 
to know how galaxies are orientated in their triaxial dark 
matter haloes. 

The orientation of galaxies in their dark matter haloes 
has received considerable attention in recent years. Studies 
of the distribu tion of sate l lite g a laxies around centrals 
in th e SDSS (|Wang et ahl 120081: lAgustsson fc Bralnerdl 
20101) and in num erical simulations (|Kang et ahl 120071 



Deason et ahl 1201 If ), studies based on the ellipticity 



corre l ation functions of galaxies (jFaltenbacher et ahl 
120091 : lOkumura et~aTI l2009f) . and studies based on angu- 
lar momentum c onsiderations in numerical simulations 
(|Bett et al.ll2010l: iHahn et al.ll2010f) all point in a similar 
direction: on average, red galaxies are aligned with their 
dark matter hosts, but with a considerable scatter between 
the position angles with a value in the range ~20-40 
degrees. The dispersion for blue galaxies is even larger. 

Scatter between the position angles irrevocably leads 
to a reduction in the anisotropy of the lensing signal; 
with lensing, we only measure the component of the dark 
matter halo ellipticity that is aligned with the lens light. 
To study the magnitude of the reduction, we use the 
simulations from the previous section, where the position 
angle of the dark matter halo is given by the position 
angle of the light, plus a Gaussian with zero mean and a 
certain width. In Figure fT^l we show the average ff°^ and 
(/ — A5) as a function of ag, the dispersion of the position 
angle distribution. We find that a oq in the range between 
20 to 40 degrees leads to a reduction in the anisotropy 
of ^25-65%. This result indicates that the scatter in 
the relative position angle significantly dilutes the shear 
anisotropy measurements. 

Qua l itative ly similar results have been obtained 
by Bettl (120111 ). who used the Millennium simulation 



(|Springel et al. 120051 ) in combination with semi-analytic 
galaxy formation models to predict the stacked projected 
axis ratio, q, of large numbers of haloes. Various alignment 
models were tested, and different methods used to measure 
the halo shapes. For most scenarios, q turned very close to 
unity, implying an nearly isotropic shear signal. 

In reality, the alignment between the position angle 
of the galaxy and the dark matter may be scale de- 
pendent. Besides that, the dark matter haloes arc not 
rigid once formed; their shapes continue to evolve in a 
way that depends on, amongst others, their formation 
history and environment; the galaxies in the centre of 
the haloes may evolve differently. Furthermore, the cen- 
tral galaxy i s expected to spherical is e the dark matter 
haloes (e.g. iKazantzidis et ahl 120101: lAbadi et all 120101: 
iMachado fc Athanassoula 2010f ). Indications already 
exist that the alignment depends on lens lumin osity 
(jFaltenbacher et ahl l2009t IMandelbaum et al.l l2006afl and 
environment ( Hahn et al. H2010T) : it could depend on other 
lens characteristics (e.g. lens ellipticity, redshift) as well. 
Finally, once more, multiple deflections and clustering 
are expected to have some impact on the anisotropy of 
the lensing signal. All these effects have to be accurately 
modeled and well understood before we can interpret any 
measured anisotropy in terms of the average property 
of dark matter haloes. Note that the study of intrinsic 
alignments of galaxies are similarly affected, and require 
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Fig. 12. Anisotropy of the lcnsing signal averaged between 
1 and 200 arcsec, for a Gaussian distributed position angle 
difference between the dark matter and the light distribu- 
tion with zero mean and width ag. Large values for ag as 
have been reported in the literature lead to significant re- 
ductions in the shear anisotropy. 



this knowledge as well for a correct interpretation. 

Our results underline the need for photometric redshifts 
- and consequently luminosities - for the lenses. Without 
photometric redshifts, we can only select lenses based on 
their colours and magnitudes. To achieve sufficient signal- 
to-noise in order to obtain competitive constraints on the 
average halo ellipticity, we have to select large numbers 
of galaxies that cover a broad range of luminosities and 
redshifts. If the average halo ellipticity depends on the 
luminosity of a galaxy (as the results from Mandclbaum 
et al. 2006a suggest), the signal-to-noise of the shear 
anisotropy measurements decrease, and may even average 
out in the worst case scenario. If luminosities are available, 
we can not only select lenses in narrow luminosity ranges, 
but also weigh the lensing measurement with luminosity, 
which improves the signal-to-noise of the lensing measure- 
ment. The lack of photometric redshifts also forces us to 
stack the lensing signal as a function of angular separation, 
rather than physical, which decreases the signal-to-noisc 
as well. This is particular disadvantageous as the shear 
anisotropy signal of an elliptical NFW profile drops very 
rapidly with increasing radius, and the signal may be 
smeared out and become undetectable. We note that a 
preliminary photometric redshift catalogue exists for the 
RCS2 for the area that has also been observed in the 
i'-band, but it only covers the redshift range z > 0.4 due 
to the absence of observations in the u-band, which limits 
its usefulness for this study. 

A new technique has recently been proposed to improve 
halo ellipticity measurements: the use of a higher order 



disto r tion of lensing known as flexion (lHawken fc Bridld 
120091: lEr fc Schneider] 120111: lEr et all 1201 lallbT l. Although 
the measurement of the flexion signal is difficult for galaxy- 
scale potentials, the first positive det ections have already 
been reported (|Velander et al.l 1201 If) . Using mock simu- 
lations o f cluste rs with SIE and elliptical NFW profiles, 
lEr et al.l (|2011bT ) find that flexion is more sensitive to the 
halo ellipticity than the shear; this may be true as well 
for stacked galaxy potentials. Furthermore, the systematic 
errors in flexion measurements differ from those in shear. 
Hence we anticipate that additional useful constraints can 
be obtained with flexion. 

6. Conclusion 

We present measurements of the anisotropy of the 
weak lensing signal around galaxies using data from the 
Red-sequence Cluster Survey 2 (RCS2). We define three 
lens samples: the 'all' sample contains all galaxies in the 
range 19 < m r / < 21.5, whereas the 'red' and 'blue' 
samples are dominated by massive low-redshift early-type 
and late-type galaxies, respectively. To study the environ- 
mental dependence of the lensing signal, we also subdivide 
each lens sample into an isolated and clustered part, and 
analyse them separately. 

We address the impact of several complications on the 
shear anisotropy measurements, including residual PSF 
systematics in the shape catalogues, multiple deflections, 
the clustering of lenses, and correlations between their 
intrinsic shapes. We run a set of idealised simulations to 
estimate the impact these might have on real data, and 
find them to be small, but not entirely negligible. We 
demonstrate that the impact of these complications can 
be reduced by a careful selection of the lens sample, i.e. 
low-redshift, massive and elliptical galaxies, as has been 
done in this work. 

We also measure the distribution of physically asso- 
ciated galaxies around the lens samples. We find that 
these satellites predominantly reside near the major axis 
of the lenses. The results of the 'red' sample are in good 
agreement with previously reported values, whilst the 
constraints of the 'all' and 'blue' sample cannot be easily 
compared as they consist of a mixture of early-type and 
late-type galaxies. 

The shear anisotropy is quantified by the anisotropy of 
the galaxy-mass cross-correlation function, (/ — f^), and 
by the ratio of the projected dark matter halo ellipticity 
and the observed galaxy ellipticity, fh- For the 'all' sample 
we find that (/ - / 45 ) = 0.23 ± 0.12, and f h = 1.50±f;gf 
for an elliptical NFW profile, which for a mean lens 
ellipticity of 0.25 corresponds to a projected halo ellipticity 
of eh — 0.38^o;25 if the halo and the lens are perfectly 
aligned. Note that various studies indicate that this may 
not be the case. These constraints provide weak support 
that galaxies are embedded in, and preferentially aligned 
with, triaxial dark matter haloes. For isolated galaxies, 
the average shear anisotropy is larger than for clustered 
galaxies; for elliptical NFW profiles, we find fh = 4.73+205 
and fh = 0.90~^{'\ 7 5 , respectively. The decrease of the 
lensing anisotropy signal around clustered galaxies may 
be due to the stripping of dark matter haloes in dense 
environments. 

We do not detect a significant shear anisotropy for 
the average 'red' lens. The shear anisotropy for the most 
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elliptical 'red' galaxies is marginally positive, but the pat- 
tern of the signal suggests that this is not the result of an 
alignment between the dark matter haloes and the galaxies. 
For the 'blue' lenses, we find that the shear is marginally 
negative, with slightly increased significance for the most 
elliptical galaxies, suggesting an anti-alignment between 
the galaxy and the dark matter. Our measurements 
highlight the need for (photometric) rcdshifts in lensing 
studies. In order to reach sufficient signal-to-noise that 
enable competitive constraints on the shear anisotropy, 
we have to stack large numbers of galaxies that span a 
broad range in luminosities and redshifts. This smears out 
the shear anisotropy, and in the worst case the anisotropy 
might even average out. 
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Appendix A: Lens selection 

We define our 'red' and 'blue' lens samples using the Wl 
photome tric redshift catalogue fr om the CFHTLenS collab- 
oration (Iffildebrandt et alJr20TlT) . This catalo gue contains 
the apparent magnitudes of galaxies observed in the same 
filters as in the RCS2, but also their photometric redshifts 
and absolute magnitudes. We split the sample in quies- 
cent (typically early-type) and star-forming (typically late- 
type) galaxies according to their photometric type identifier 
Tbpz- 2bpz corresponds to the galaxy type of the best- fit 
template (E/S0, Sbc, Scd, Im, SB3, SB2, with T BPZ values 
from 1.0 to 6.0). We select the early-types and late-types 
with Tbpz < 1-5 and Tbpz > 1-5 respectively. These selec- 
tions are further divided into a bright and faint sample us- 
ing M r i < —22.5 and M r i > —22.5. For these four samples, 
we plot the colours as a function of magnitude in Figure 
IA.1I Wc define selection boxes of the 'red' and 'blue' lens 



sample, which are aimed at selecting the bright early-types 
and late-types. Figure IA~T1 shows that the 'red' lens sample 
contains almost no late-type galaxies, nor faint early-types. 
Unfortunately, for the 'blue' sample we cannot exclusively 
select luminous late-type galaxies, and the sample also con- 
tains faint late- type and early-type galaxies. These selection 
criteria are applied to the catalogues of the RCS2 to select 
the 'red' and 'blue' lens sample. 

To study how well we can separate early-types from 
late- types, we compare our selection to previously employed 
separation criteria. We find that 98% of the 'red' lenses, 
57% of the 'blue' lenses and 28% of the 'all' lenses have a 
photometric type Tbpz < 1-5, and therefore have a spec- 
tral energy distributions similar to red elliptical galaxies. 
The 'red' sample therefore barely contains late-type galax- 
ies. The majority of the 'blue' sample are actually faint red 
early-type galaxies, as can be seen from Figure IA.11 The 
lensing signal is dominated, however, by the massive late- 
type galaxies. The faint lenses mainly add noise. The purity 
of the 'blue' sample could be improved by shifting the se- 
lection boxes to bluer colours, but this at the expense of 
removing the majority of massive late-type lenses. Finally, 
the majority of the 'all' sample are blue and not very mas- 
sive late-type galaxies. As a check, we also compare to the 
u — r colour select i on crit erion, which has been used in 
iMandelbaum et al~l (|2006al ). In this work, galaxies with an 
SDSS u — r > 2.22 model colour are selected for the red 
sample, whilst galaxies with u — r < 2.22 are selected as 
blue galaxies. When we select galaxies based on identical 
u' — r' criteria (hence ignoring small differences between the 
filters), we find very similar results: 99% of the 'red' lenses, 
58% of the 'blue' lenses and 30% of the 'all' lenses have a 
colour u' — r' > 2.22 and are red. 

Appendix B: Environment selection 

We subdivide the lenses in a clustered and an isolated sam- 
ple, depending on whether or not the lens has a neighbour 
within a certain projected radius range that has a lower 
apparent magnitude than the lens. To determine which 
radius effectively separates the lenses into a low-density 
and a high-density sample, we compare the lensing signals 
for four lens selections: those that have a brighter neigh- 
bour within 30 arcsecs, those with a brighter neighbour 
between 30 arcsecs and 1 arcmin, those with a brighter 
neighbour between 1 arcmin and 2 arcmins, and those that 
are the brightest object within 2 arcmins. We show the 
lensing signal of the four selections in Figure IB.ll We find 
that the lensing signal at scales smaller than ~0.5 arcmin 
does not change much using the different isolation criteria, 
which demonstrates that we select haloes of similar mass. 
However, the large-scale signal decreases significantly for 
an increasing minimum separation to a brighter neighbour. 
The galaxies with no brighter object within 1 arcmin are 
selected for the isolated sample, the other galaxies are 
selected for the clustered sample. We could in principle 
select a more clearly distinguished sample of isolated and 
clustered lenses, e.g. by only selecting those galaxies with 
no brighter neighbour within 2 arcmins, and those with a 
brighter neighbour within 30 arcsecs. However, this would 
reduce the signal-to-noise of the lensing measurements 
such that no useful constraints could be obtained on the 
average halo ellipticity, and we therefore choose not to. 
Note that the dependence of the large-scale lensing sig- 
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Fig. A.l. Colour- magnitude diagrams that are used to define the lens sample selection. In gray we show the loca- 
tions of the galaxies wit h m r > < 21.5 from the Wl photometric redshift catalogue from the CFHTLenS collaboration 
(jHildebrandt et al.l[201lt) . In the top row we show the late-type galaxies with an absolute magnitude M r i < —22.5, in 
the second row we show the late- types with M r > > —22.5, and in the third and fourth row we show the early- types 
with M r i < —22.5 and M r > > —22.5, respectively. The coloured boxes illustrate the lens sample selection criteria; all 
galaxies that reside in all the red dashed (green solid) boxes form the 'red' ('blue') lens sample. We find that the 'red' 
lens sample contains very few late- types, as they are excluded in the g' — r 1 versus g 1 diagram (first column). Also, most 
faint early-type galaxies are excluded from this sample as well (see, e.g., the third and fourth columns). For the 'blue' 
lens sample, we cannot define selection criteria that exclusively select bright late-type galaxies, and this sample therefore 
also contains a number of faint late-type and early-type galaxies. 



nal on the distance to a brighter neighbour is partly caused 
by differences in the lens environment, and partly by 
differences of the projected densities along the line-of-sight 
(LOS). For the clustered lenses, the LOS projections do 
not simply add random noise, but systematically increase 
the large-scale lensing signal because of the neighbour se- 
lection. Similarly, for the isolated galaxies, the underdense 
LOS leads to a decrease of the lensing signal on large 
scales. For the purpose of this work this is not important, 
since our main focus is to measure the anisotropy of the 
lensing signal. However, it should be kept in mind in the 
interpretation of the results. 



Appendix C: Average ratios and their errors 

To calculate / mm , Z,™^ and (/ — f^) we have to determine 
the ratio of two variables. Let us call these variables x and 
y, and m = x/y the ratio we are interested in. If x and 
y are independent and have a Gaussian distribution with 
means fx x and [i y and widths a x and g„ we can com pute 
the probability distribution of m using (lHinklevlll96"9h 
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Fig. B.l. The lensing signal as a function of projected separation for the three lens samples that have been divided 
according to the minimum distance to a brighter galaxy, as indicated in the top right panel. The red diamonds indicate 
the lensing signal of those lenses that have a brighter neighbour within 30 arcsecs, the green squares the signal of the 
lenses that have a brighter neighbour between 30 and 60 arcsecs, the blue triangles the signal of the lenses that have 
a brighter neighbour between 60 and 120 arcsecs, and the orange circles the signal of the lenses that have no brighter 
neighbour within 2 arcmins. The dashed lines show the best fit NFW profiles fitted to the lensing signal of the total 
sample between 50 and 500 kpc using the mean lens redshift. The lensing signal on scales <0.5 arcmin is roughly similar 
for all lens samples, but the lensing signal on larger scales clearly decreases for increasingly isolated galaxies. 
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By integrating p(m) we can determine the median, and the 
68% confidence intervals. 

Next, we want to combine various measures of rrij 
into one average. We fol low the approach described in 
iMandelbaum et al . (2006a). It relics on the use of a slightly 
different variable, i.e. yi — mxi, which is again a random 
Gaussian variable. If the shear ratio is constant over the 
range of interest, its value can be determined with 



JywiiVi - mx i) < 
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id, 



(C.6) 



where 



and <jy, and at. the error on 



Wi = l/(cr^. + m 2 c 
yi and Xi, respectively. Z = then gives the average ratio, 
whilst Z = 1 gives the 68% confidence intervals. Note that 
Equation (|C.6I) can also be used to determine the ratio of 
a single measurement, but we prefer the use of Equation 
(|C.1|) as it directly provides the full probability distribution. 



Appendix D: Lens light contamination 

The light from bright and elliptical lenses changes the 
source number density along the major and minor axis dif- 



ferently on small projected separations close to the lens. 
This could bias the correction we make to account for phys- 
ically associated galaxies in the source sample. We inves- 
tigate the size of the effect by selecting all galaxies with 
16 < m r / < 18, 18 < m r > < 19 and 19 < m r ' < 20, and di- 
vide each selection in three ellipticity bins. For each of these 
samples, we measure the source number density in the ma- 
jor and minor axis quadrants, and show the results in Figure 
ID. 11 On small scales, we find that the deficiency in the ma- 
jor axis quadrants is indeed significantly larger, and the 
difference increases for brighter and more elliptical galax- 
ies. For projected separations larger than 0.2 arcmin, the 
difference is negligible, except for the brightest and most 
elliptical bins. Therefore, for galaxies with m r < < 19 we 
only use scales larger than 0.2 arcmin, whilst for galaxies 
with m r > > 19 we use scales larger than 0.1. 



Appendix E: Magnification 

Magnification may also affect the shear anisotropy measure- 
ments. In this appendix, we address two potential contami- 
nants: anisotropic magnification by the lens, and additional 
magnification by other foreground galaxies. 



E.l. Anisotropic magnification 

If dark matter haloes are triaxial, the magnification of back- 
ground sources will have an azimuthal dependence. Hence if 
the galaxy and dark matter halo arc aligned, part of the ob- 
served anisotropy in the distribution of source galaxies may 
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Fig. D.l. The source number density along the major and 
minor axis close to bright foreground galaxies, indicated by 
green squares and blue triangles, respectively, selected on 
their apparent magnitudes (columns) and ellipticity (rows). 
On small scales, the deficiency of background galaxies along 
the major axis is larger for brighter and more elliptical 
galaxies. 



be due to anisotropic magnification. The shear anisotropy 
measurements are corrected, however, assuming that the 
anisotropic distribution is solely due to physically associ- 
ated galaxies. If anisotropic magnification has a strong ef- 
fect on the source density, the correction might therefore be 
biased. To obtain an estimate of the impact of anisotropic 
magnification, we assume that the galaxy and the dark mat- 
ter halo are perfectly aligned. Furthermore, we assume that 
at small scales the stacked density profile of the lenses is 
approxim ately described by an SIE, with a surface density 
given by (jMandelbaum et al"1l2006ah 



W DjDu 1 
c 2 £>, 2r 



1 + y cos(2(9) 



(E.l) 



a is the average velocity dispersion, which we determine 
by fitting an SIS to the azimuthally averaged tangential 
shear within the virial radius, Di s /D s is determined by in- 
tegrating over the source redshift distribution, and the halo 
ellipticity is assumed to be equal to the mean galaxy el- 
lipticity, as tabulated in Table [1] (i.e. assuming f\ x = 1). 
We determine the magnification using fx = 1 + 2k at a 
projected separation of 10 arcsecs. The change in number 
density due to magnificatio n in the major axis quadrant, 
/^ as (r), is calculated with (lNaravanlll989l) 
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where No(m) is the background galaxy number density, 
N-Q(m,r,6) is the number density after magnification in 
the major axis quadrant, and s(m) is the slope of the 
logarithmic galaxy number counts at magnitude m. The 
change in number density along the minor axis, /™ ag (f), 
is calculated similarly. We det ermine s(m) usin g the 
photometric redshift catalogue of lllbert et all (|2006t ) , and 
find that it decreases from 0.34 at m r > = 22 to 0.22 at 
m r ' = 24; we use the average value s(m) — 0.28. Hence 
the number density of our source sample is diluted due 
to magnification, more strongly along the major axis of 
the lenses than along the minor axis as the magnification 
along the major axis is larger. For the 'all' lens sample, we 
find that at 10 arcsecs, /" ag = 0.9969 and /™ ag = 0.9974. 
For the 'red' lens sample we find /g lag = 0.9927 and 
/!f ag = 0.9935, and for the 'blue' lens sample we find 
/™ a s = 0.9946 and /™ ag = 0.9954. Hence the effect of 
anisotropic magnification is expected to be very small, 
and does not cause the observed anisotropy of the excess 
source galaxy density ratio. 



E.2. Magnification by L2 

Foreground galaxies that lens both the lens galaxy 
and the sources also magnify the background sky. This 
causes a change in the source number density around the 
foreground galaxy. Consequently, the positions of the fore- 
ground galaxy and the background sources become related. 
This could potentially lead to a false shear anisotropy 
signal. We explain the effect using a cartoon in Figure [ETl 
In this figure, we measure the shear anisotropy around 
lens LI. The intrinsic shape of LI is shown by the black 
ellipse, and the regions where we average the tangential 
shear in the absence of a foreground lens are indicated by 
the black lines. In the presence of a foreground lens L2 
that is located close to one of the original quadrant's axes, 
the position angle of LI changes, and the quadrants where 
we average the shear rotate, as indicated by the red dotted 
lines. L2 also locally magnifies the background sky, leading 
to either an increase or decrease in the source number 
density around L2. When we average the tangential shear 
in the rotated quadrants as would be done in observations, 
these magnified sources predominantly move into the 
minor axis quadrant (region A in Figure [3]) . If the source 
number density increases around L2, we find that we 
have more sources in the minor axis quadrants than in 
the major axis quadrants. Furthermore, the average shear 
changes, biasing the shear ratio high. If the source number 
density decreases instead, the opposite effect happens. 

The presence of this bias could be identified by 
comparing the number counts of the sources in the two 
quadrants for the radial bins close to the lens, as the 
effect is strongest on small scales. The effect is mixed, 
however, with the anisotropic distribution of satellite 
galaxies. However, we already found in Appendix IE. II 
that magnification for our source sample is negligible. 
Therefore, we expect that this source of bias is small, and 
can be ignored with the current data. The effect may be 
measurable by selecting a sample of source galaxies that 
are clearly in the background, and whose number density 
slope is steep (e.g. Lyman-break galaxies), as these are the 
conditions favourable to magnification. 
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Fig. F.l. Observed shear ratio of an elliptical lens LI at z = 0.4 for sources at a projected separation of 100 arcsecs, 
in the presence of an additional foreground galaxy L2 at z = 0.1 that is located at (x,y). In the left-hand panel we 
show /mm; the ratio of the average shear along the major and minor axis. The location of LI is indicated by the blue 
ellipse in the centre of the image, the location of the sources by the cyan dots. We find that depending on the position 
of L2, the shear ratio either increases or decreases. The black solid lines indicate the position of L2 where / mm remains 
constant, the red dot-dashed (green dashed) lines indicate a decrease (increase) of / mm of 0.2, 0.4 and 0.6, respectively. 
In the right-hand panel, we show /™mi the shear ratio corrected using the cross terms. For large L2-source separations, 
the induced shear on the lens and the sources is almost constant, and is effectively removed using the cross terms. For 
small separations, the residual is large as the applied shear from L2 varies strongly along the ring of sources. Note that 
both LI and L2 are modeled by an SIS. 



Appendix F: Multiple deflections 

To visualize how the shear anisotropy is affected by multiple 
deflections, we simulate a lens galaxy LI at a redshift 0.4. 
We assume the LI is an SIE, with an ellipticity ei = 0.2 
and e2 = 0. The dark matter halo is perfectly aligned with 
the light distribution, and has the same ellipticity as the 
lens. We compute the lensing signal with Equation (|2~i|) 
for a velocity dispersion of 200 km s _1 . In the absence of 
multiple deflection, we find that / mm = (it + e/j) / (it — eh) 
1 . 1 .'iti. We insert a second lens L2 in the image, a round 
SIS at z = 0.1 with a velocity dispersion a = 200 km 
s _1 . For every position of L2 in the simulated image, we 
calculate / mm (r = 100), the shear ratio we would observe 
for a projected lens-source separation of 100 arcsecs, using 
a source redshift of 0.8. The result is shown in Figure lF4l 

We find that depending on the location of L2, / mm either 
becomes larger or smaller than 1.136, and even in some 
configurations turns negative. These trends can easily be 
understood: e.g. if L2 is located at (50,250), left to the lens 
and the ring of sources, it increases the tangential shear 
of the sources along the major axis (region B in Figure [3]), 
but it decreases the tangential shear of the sources along 
the minor axis (region A in Figure [3]). The impact on the 
ellipticity of LI is very small, and consequently we find that 
/mm increases. If L2 is located very close to the sources, it 
can change the source ellipticities by such amounts that the 
net tangential shear with respect to LI becomes negative, 
which results in a negative shear ratio. 

Equivalently to Equation (|23[) . we compute as a 



function of the position of L2, which is shown in the right- 
hand panel of Figure IF. 1 1 If L2 is at a large distance from 
the sources and LI, the induced shear on both LI and the 
sources is almost constant, and is effectively removed using 
the cross terms. For smaller separations, correcting the ratio 
using the cross terms does not work well as the applied 
shear from L2 varies strongly along the ring of sources. 

To quantify the net effect of multiple deflections on the 
shear ratio, we determine the total contribution to "ft, a and 
It.B by integrating over all L2 positions. Then we compute 
the average value of the ratio: 

(/rnrnM) = ^4^ , (F.l) 

7t,A0 + n L2 i\"ft,A 

where "ft,Ao and "ft, bo are the average tangential shear along 
the minor and major axis in the absence of L2, A.j t ,A and 
Ajt.B are the total contributions to the tangential shear 
along the minor and major axis, and n~L2 is the foreground 
galaxy number density. The change of the shear ratio de- 
pends on the ellipticity of the lens, the density profiles of 
LI and L2, the area within which we integrate the contri- 
butions of L2, and on the number of second deflectors, n~i,i. 
We study these dependencies below. 

In Figure lFT2k . we show (/ mm ) and (/^ r ) as a function 
of r ou t, the radius of the circle centred at LI within which 
we integrate the contributions of L2 (hence for a r out of 50 
arcsecs, we only account for contributions of L2 that are 
located within 50 arcsecs from LI), for a Ll-sourcc separa- 
tion of 100 arcsec. We use three different profiles for L2, i.e. 
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Fig. F.2. (a) Observed value of / mm (dotted line) and /£J^ (solid line) for a lens LI with an SIE profile and ellipticity 
(ei,e2) = (0.2,0.0), measured using sources at 100 arcscc, as a function of r out , the projected radius from LI within 
which the contributions of L2 are integrated. The black lines show the reduction if L2 is an SIS, the green (red) if L2 is 
a TIS with a truncation radius of 150 (50) arcsec. The dashed line indicates the shear ratio in the absence of L2. This 
figure demonstrates that L2 at projected separations > 100 arcsec from LI have a non-negligible influence on the shear 
ratio. Multiple deflections reduce the observed shear ratio, but using the cross shear we can correct for most of it. (b) 
Observed value of f™™ as a function of source radius, for an r out of 4000 arcsec. The reduction of f^™ increases with 
increasing Ll-source separation, (c) Difference between the observed value of f™m and the shear ratio in the absence of 
multiple deflections (/ mm . n o md) as a function of the ellipticity of LI, measured at a Ll-source separation of 100 arcsec, 
for an outer radius of integration of 4000 arcsec. The reduction of the shear ratio increases with decreasing ellipticity. 




Fig. F.3. Net effect on /A£j so (left-hand panel) and (/ — /45)AE; so (right-hand panel) for a LI SIS with an ellipticity 
(ei,e2) = (0.2,0.0), measured using sources at 100 arcsec, in the presence of a second foreground galaxy L2 located at 
(x, y). The location of LI is indicated by the blue ellipse in the centre of the image, the location of the sources by the 
cyan dots. We find that depending on the position of L2, the contribution to the shear is either positive or negative. The 
black solid lines indicate the position of L2 with no net contributions to /A£i so and (/ — /45)AE iso , the red dot-dashed 
(green dashed) lines indicate a decrease (increase) of AE of 1, 2 and 3 [M /pc 2 ] respectively. 



an SIS and two truncated SIS (TIS) profiles with a trunca- 
tion radius of 50 and 150 arcsec. Most real galaxies have a 
density distribution that falls somewhere in between these 
extremes. The value for n~t2 we have adopted is 100 galaxies 
per square degree. The velocity dispersions and redshifts of 
the galaxies are similar to the values used before. 



We find that the reduction of the shear ratio is a strong 
function of r ou t- It reaches a minimum when r out is equal 
to the Ll-source projected separation, which is expected 
from Figure [FT] if L2 is located at the ring of sources, the 
shear ratio becomes large and negative, which leads to a 
reduction of the average ratio. For larger r out , we find that 



26 



Edo van Uitert et al.: Constraints on the shapes of dark matter haloes 




Sources 



Appendix G: Intrinsic alignments 

The shear anisotropy (/— /45)AEi SO of the 'red' lens sample 
on small scales is negative (Figure [TU]) . Since the source 
sample contamination of physically associated galaxies is 
largest at the same scales (see Figure 0|, and since the 
distribution of these satellite galaxies is anisotropic (see 
Figure [5]), a negative shear signal could also be caused if 
these satellite galaxies are preferentially radially aligned. 
To estimate the value of the average tangential intrinsic 
alignment that would produce such a signal, we rewrite 
Equation ©: 



.4 



/A£ iso (r) = — J2 Wi&EieZi cos(2A^)- 



N 



(G.l) 



5^ic i AX3 orit 7 / e«iCos(2Afl i ); 



N 



4=1 



wiefc cos 2 (2A^) 



N, 



2Yw l e 2 a l cos 2 (2A6 i ), 



Fig. E.l. Cartoon to illustrate a potential bias in the halo 
ellipticity measurement due to magnification. In the pres- 
ence of a foreground galaxy L2, the quadrants where the 
tangential shear is averaged are rotated. L2 also magnifies 
the background sky, leading to an relative increase or de- 
crease in number of sources in the minor axis quadrant. 
This potentially biases the measured shear ratio on small 
scales. 



the shear ratio increases again. If L2 is an SIS, we find that 
the ratio turns over and continues to decrease. If L2 is a 
TIS, it does not contribute to the shear ratio if it is located 
at a projected separation much larger than the truncation 
radius, and the ratio therefore converges to a certain value. 
We observe that the impact of multiple deflections is mostly 
removed for (fmm)- Even if L2 is an SIS, the total reduc- 
tion of the shear ratio is small as long as we integrate the 
contributions of L2 over a sufficiently large area. 

In Figure lR2b . we show (/mm) and (f™m) as a function 
of projected source separation for an r out of 4000 arcsec. We 
find that (/mm) decreases with r s if L2 is an SIS. If L2 is 
a TIS, (/mm) reaches a minimum because L2 cannot shear 
both LI and the sources if r s becomes much larger than 
the truncation radius, which reduces the impact on (/ m m)- 
Furthermore, we find that (/m°m) is hardly affected by mul- 
tiple deflections on small scales. On large scales, however, 
the contribution from L2 is no longer constant, and the cor- 
rection scheme fails. 

Finally, we show the bias of the shear ratio as a function 
of ellipticity of LI in Figure lFT2b . We find that the reduction 
is largest for the roundest lenses, as their position angles 
are affected most by the presence of L2. The impact of mul- 
tiple deflections decreases for lenses with larger cllipticities. 

For completeness, we show the impact of an addi- 
tional foreground lens L2 located at (x, y) on /AS iso and 
(/ — ,/45)AS iso in Figure IF. 31 Similarly as for the shear 
ratios shown in Figure IF. II we find that the cross terms 
remove most of the systematic contributions as long as L2 
is not located very nearby. 



with 7/ the tangential intrinsic alignment, and N s and 
Ni the number of sources and physically associated galax- 
ies, respectively. We assume that the shear signal is 
isotropic, hence the first term on the right-hand side of 
Equation (|G.1[) cancels. As the shear anisotropy is ex- 
pected to be positive for the red early-type galaxies that 
make up the 'red' lens sample, our estimate of 7/ is a 
lower limit. To account for the fraction of satellites that 
is anisotropically distributed, we multiply Equation (jG.ll) 
with f cg (r, A#)/(/ cg (r))AG- For / cg we adopt the same form 
as in Section \%~M i.e. Equation (|15|) . which we insert into 
Equation 



(jG.lll . Ignoring the isotropic part, which cancels 
when averaged over the angle, we obtain: 

A Nl IN 
/AE iso (r) = - ]T m i AS crit7 /e^ cos 2 (2A0 i ) ' 



and therefore 
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/A£ iso (r) N s + Ni N iso - 1 
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Ni N Ae 
/AE iso (r) N is 
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(G.2) 
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2AT, 



crit 



A Ae ' 



where we used that (N s + A/)/A/ = A iso /(A iso - 1). 
Averaging the results of the first two radial bins of the 
'red' lenses, we obtain 7/ = — 0.062l ; 062 for a = 1 at an 
average separation of 0.19 arcmin (~65 kpc at the mean 
lens rcdshift). Our constraints are not particularly com- 
petitive (compare, e.g., Hirata et al. 2004), mainly because 
our methods are not designed for this measurement. For 
example, had we included brighter galaxies in the source 
sample, a larger fraction of those would be physically 
associated to the lenses, improving the constraints on the 
galaxy overdensity, and hence on 7/. This, however, is not 
the purpose of this work. Note that the value and error 
on 77 we o btain is roughly a factor 10 larger than the 
results from iHirata et al.l (|2004f ) , who aimed their analysis 
to measure this effect using data from the SDSS (although 
no division was made between early-type and late-type 
galaxies, which could average out the effect). Conversely, 
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the expected impact of intrinsic alignments is about a 
factor ten smaller than the signal we observe, hence it is 
unlikely that our measurements are significantly affected. 
Hence the negative shear anisotropy is unlikely caused by 
intrinsic alignments. 
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